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ABSTRACT 

O \ We present high-resolution (FWHM ~ 0.3-1.5 km s~^) spectra, obtained with the 

\ AAT UHRF, the McDonald Observatory 2.7-m coude spectrograph, and/or the KPNO 

' coude feed, of interstellar Ca I absorption toward 30 Galactic stars. Comparisons of 

the column densities of Ca I, Ca II, K I, and other species — for individual components 
^ ■ identified in the line profiles and also when integrated over entire lines of sight — yield 

. information on relative electron densities and depletions (dependent on assumptions re- 

garding the ionization equilibrium). There is no obvious relationship between the ratio 
iV(Ca l)/A^(Ca II) [equal to ne/(r/ar) for photoionization equilibrium] and the fraction 
of hydrogen in molecular form /(H2) (often taken to be indicative of the local density 
wh)- For a smaller sample of sightlines for which the thermal pressure {n-^T) and local 
density can be estimated via analysis of the C I fine-structure excitation, the average 
electron density inferred from C, Na, and K (assuming photoionization equilibrium) 
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seems to be independent of tih and heT. While the electron density (ng) obtained from 
the ratio A'^(Ca l)/A^(Ca II) is often significantly higher than the values derived from 
other elements, the patterns of relative rig derived from different elements show both 
similarities and differences for different lines of sight — suggesting that additional pro- 
cesses besides photoionization and radiative recombination commonly and significantly 
affect the ionization balance of heavy elements in diffuse interstellar clouds. Such addi- 
tional processes may also contribute to the (apparently) larger than expected fractional 
ionizations (ne/nn) found for some lines of sight with independent determinations of 
njj- In general, inclusion of "grain- assisted" recombination does reduce the inferred n^,, 
but it does not reconcile the ng estimated from different elements; it may, however, 
suggest some dependence of rig on nn- The depletion of calcium may have a much 
weaker dependence on density than was suggested by earlier comparisons with CH and 
CN. Two appendices present similar high-resolution spectra of Fc I for a few stars and 
give a compilation of column density data for Ca I, Ca II, Fe I, and S I. 

Subject headings: ISM: abundances, ISM: atoms, ISM: kinematics and dynamics. Line: 
profiles 

1. Introduction 

High-resolution (FWHM ~ 0.3-1.8 km s~^) surveys of the optical interstellar absorption lines 
of Na I (Welty, Hobbs, & Kulkarni 1994; hereafter WHK), Ca II (Welty, Morton, & Hobbs 1996; 
hereafter WMH), and K I (Welty & Hobbs 2001; hereafter WH) have revealed complex velocity 
structure in many Galactic lines of sight. Fits to the absorption-line profiles have yielded column 
densities, line widths, and velocities for several hundred individual components in each case — 
enabling statistical studies of the properties of the corresponding individual interstellar clouds. 
Under the usual assumption of photoionization equilibrium, comparisons of such high-resolution 
spectra of K I and Ca II with similar spectra of Ca I can provide information on the relative 
electron densities and elemental depletions in the individual components. This information can 
then be used to model the profiles of many other neutral and singly ionized species seen in lower 
resolution UV spectra — yielding detailed abundances, depletions, and physical conditions for the 
individual clouds (e.g., Welty et al. 1999b). 

Because even the strongest line of Ca I (at 4226. 7A) is typically quite weak (W^ < 10 mA), 
however, very few spectra at resolutions better than 3 km s~^ have ever been reported (Hobbs 1971; 
Welty et al. 1999b). The earliest quantitative comparisons of Ca I and Ca II yielded mean electron 
densities (rig) ranging from about 0.1-1.3 cm~^ (Hobbs 1971; White 1973; Federman &; Hobbs 
1983; but using the slightly different atomic data adopted in this study). If a fractional ionization 
ng/nfj comparable to the interstellar carbon abundance is assumed (i.e., if photoionization of C I is 
the primary source of electrons) , such values of rig would imply fairly high local hydrogen densities 
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(nn) and thermal pressures (nnT). Alternatively, if the true nn are lower, then those Ue might 
suggest that hydrogen is slightly ionized even in predominantly neutral clouds. 

While ratios such as N{Ca I)/N{Ca II) can in principle yield estimates for rtg, recent observa- 
tional and theoretical studies have suggested that additional processes besides photoionization and 
radiative recombination may commonly and significantly affect the ionization balance of heavy ele- 
ments in diffuse clouds. For example, the electron densities derived from different A^(X I)/A^(X II) 
ratios in a given line of sight can differ by factors of 10 or more (Fitzpatrick Sz Spitzer 1997; 
Welty et al. 1999a, 1999b). Moreover, the corresponding fractional ionizations Ue/nu (with nn 
estimated from the C I fine-structure excitation equilibrium) can be significantly higher than the 
value "expected" from photoionization of carbon. While some additional electrons may arise from 
cosmic ray and/or x-ray ionization of hydrogen (e.g., Wolfire et al. 1995), it has been suspected 
for some time that electron densities derived under the assumption of photoionization equilibrium 
may be spuriously high (e.g.. Snow 1975). Another potentially relevant process is charge exchange 
between the various dominant singly ionized species and neutral or negatively charged small grains 
(or large molecules) , which would enhance the corresponding trace neutral species — thus yielding 
spuriously high rig if only radiative recombination is assumed (Lepp et al. 1988; Weingartner & 
Draine 2001; Liszt 2003). Accurate abundances of both neutral and singly ionized species — for a 
number of elements and for sightlines sampling a variety of interstellar environments — are needed 
to identify and disentangle the various processes that may affect the ionization balance. 

In this paper, we discuss spectra of interstellar Ca I, observed at resolutions of 0.3-1.5 km s~^, 
toward 30 relatively bright Galactic stars. The high spectral resolution enables more detailed, 
component-by-component comparisons with similar spectra of Na I, K I, Ca II, and various molec- 
ular species — and thus more accurate characterization of the properties of individual interstellar 
clouds. In Sections 2 and 3, we describe the procedures used to obtain, reduce, and analyze the 
spectra. In Section 4, we discuss the statistical properties of the ensemble of individual components 
detected in Ca I, we compare the Ca I absorption with that due to other species, and we explore 
possible implications for determinations of the electron density and depletions in interstellar clouds. 
Comparisons of various A''(X 1)/N{X II) ratios for five lines of sight reveal both similarities and 
differences in the patterns of relative n,. inferred from the different elements — suggesting that 
several different processes may affect the ionization balance. In Section 5, we give a brief summary 
of the results. In two appendices, we show similar high-resolution spectra of Fe I toward six stars 
and we list available column density measurements for Ca I, Ca II, S I, and Fe I — supplementing 
the values for Htot (= H I 2 H2), H2, Na I, K I, Li I, C I, and CH tabulated by WH. 
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2. Observations and Data Reduction 

2.1. Program Stars 

The 30 stars observed in the present high-resolution Ca I survey are Usted in Table 1. All 
but one of the stars (a Cam) were included in at least one of our previous surveys (WHK, WMH, 
WH), where we have tabulated Galactic coordinates, V, E{B — V)^ spectral types, and distances. 
Because the Ca I absorption is typically quite weak, most of the stars chosen for observation in this 
survey have E{B -V) > 0.05. 

2.2. Instrumental Setups 

Because the instrumental setups employed for these observations of Ca I are very similar to 
those used for the Na I, Ca II, and K I surveys, we will give only brief summaries here; more 
extensive descriptions may be found in WHK and WMH. 

High-resolution Ca I spectra, with Av = 0.53ib0.05 km (FWHM), were obtained for 10 
stars with the McDonald Observatory (McD) 2.7-m telescope and coude spectrograph in 1996 
September. As for the Na I and K I spectra reported by WHK and WH, the echelle spectrograph 
was used in the double-pass configuration. A slit width of lOOyum, corresponding to 0.23 arcscc on 
the sky, was used for all observations. The TK-4 CCD (1024 x 1024, 24 /xm pixels) was placed at 
the "scanner" focus. The dispersion achieved with this setup corresponds to 3.1 mA per pixel, or 
2.4 pixels per resolution element. A quartz lamp was used for fiat-field exposures; a Th-Ar hollow 
cathode lamp was Tiscd to determine the wavelength scale. The resolution was determined from 
the measured widths of the thorium lines in the Th-Ar exposures, assuming an intrinsic (thermally 
broadened) width of 0.55 km s""*^ (Willmarth 1994, priv. comm.; Crawford et al. 1994; WHK). A 
signal-to-noise ratio (S/N) ~ 165 (per half resolution element) was achieved in 120 minutes (two 
60-minute exposures) for ^ Per, which has B ~ 4.1. 

Higher resolution Ca I spectra, with = 0.32 km s~^ (FWHM), were obtained for 5 stars 
with the Ultra-High Resolution Facility (UHRF; Diego et al. 1995) on the 3.9-m Anglo- Australian 
Telescope in 1997 March. The 35-slice image sheer, which takes light from a 1.0 x 1.5 arcsec 
aperture to form a 30/im wide pseudo-slit (Diego 1993). and the Tektronix CCD (1024 x 1024, 
24 ^m. pixels) were used for all the observations reported here. The dispersion achieved with this 
setup corresponds to 2.2 mA per pixel, or 2 pixels per resolution element. Because of the large 
number of CCD rows filled by the stellar spectrum and the clean response of the CCD, flat fielding 
was not required (for the S/Ns achieved in this study). The spectral resolution, A-y = 0.32 km s~^, 
was determined via scans of a stabilized He-Ne laser line at 6328.160 A and by examination of 
the widths of the thorium lines in the Th-Ar exposures used for wavelength calibration. Individual 
exposures were limited to 30 minutes to preserve the high spectral resolution; S/N ^ 225 was 
achieved in three exposures totalling 75 minutes for 5 Sco, which has B ~ 2.2. 
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Cal spectra with a resolution Av = 1.3-1.5 km (FWHM) were obtained for 22 stars with 
the 0.9-m coude feed (CF) telescope and coude spectrograph of the Kitt Peak National Observatory 
during various runs from 1993 to 2000. To achieve the desired high resolution, we used the cchelle 
grating, camera 6, and a 90-120 //m entrance slit, corresponding to 0.6-0.8 arcsec on the sky. 
Several different grism, wedge, and filter combinations were used in the various runs to cross- 
disperse the spectra and to isolate the desired spectral orders. Two different CCD's were used 
(TIKA: 1024 x 1024, 24 //m pixels in 1993; F3KB: 1024 x 3072, 15 /xm pixels thereafter). The 
dispersion achieved with this setup corresponds to 8.3 mA per pixel in 1993 and 5.2 mA per pixel 
thereafter, corresponding to about 2.4 and 3.8 pixels per resolution clement, respectively; the post- 
1993 spectra were re-binned to approximately optimal sampling. Flat-field exposures were obtained 
using a quartz lamp; the wavelength calibration was established via a Th-Ar lamp. We achieved 
S/N ~ 220 in 100 minutes for ^ Per. 



2.3. Data Reduction 

The initial processing of the CCD images from all three telescopes followed procedures similar 

to those used in our previous surveys (WHK, WMH, WH). Various routines within IRAF were 
employed to subtract the bias and divide by a normalized flat fleld (except for the UHRF data) ; the 
one-dimensional spectra then were extracted using the apcxtract package, with variance weighting. 
Wavelength calibration was accomplished via exposures of Th-Ar hollow cathode lamps, using the 
thorium rest wavelengths of Palmer & Engelman (1983). Multiple spectra for a given target were 
sine interpolated to a common heliocentric wavelength grid and summed. 

Comparisons between the individual Th-Ar exposures and between the Ca I spectra and pre- 
viously obtained high-resolution spectra of Na I, K I, and Ca II yield complementary information 

on the accuracy of the velocities. For the coude feed spectra, quadratic fits to the positions of the 
13-28 lines found in the Ca I order yield typical residuals of 0.8 mA, or 0.06 km s^^. Calibration 
exposures were usually obtained at the beginning, middle, and end of each night; in general, the 
zero points agree to within 5 mA, or 0.35 km s~^ (and usually much better), and the dispersions 
agree to within 0.02% for any individual night. At the higher McDonald and UHRF resolutions, 
only 4-5 lines with previously known wavelengths were generally found within the ~ 1.5-2.3 A 
spectral range recorded on the CCDs. Linear fits to the positions of those lines yield typical resid- 
uals of 0.3-0.8 mA, or 0.02-0.06 km s~^. For the UHRF observations, calibration spectra were 
obtained several times during each night. The zero points agree to within 5 mA, or 0.35 km s^^, 
and the dispersions agree to better than 0.05% for each night. Comparisons of the interstellar line 
profiles for seven lines of sight observed at both McDonald and the coude feed suggest that any 
systematic differences in velocity zero point are less than 0.2 km s~^. There also appears to be 
good agreement in both velocity scale and zero point between these Ca I spectra and the Na I, 
Ca II, and K I spectra reported by WHK, WMH, and WH. For example, for 35 "corresponding" 
components identified in both Ca I and K I in the profile analysis (see below), f (Ca I) — v{K I) is 
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0.02 lb 0.19 km s ^ (mean ± standard deviation). 

The summed, wavelength-calibrated spectra were normalized by fitting Legendre polynomials 
to regions free of interstellar absorption. In all cases, sufficient continuum was present on both sides 
of the interstellar features. The S/Ns determined in the continuum fits and the total equivalent 
widths of the interstellar Ca I absorption measured from the normalized spectra are given in 
columns 2-3 (CF), 4-5 (McD), and 6-7 (UHRF) of Table 1. The median S/Ns achieved for the 
three data sets are ~ 120, 175, and 215, respectively. Those median S/Ns yield 2(7 equivalent 
width uncertainties or limits for narrow, unresolved absorption lines (including contributions from 
both photon noise and continuum placement uncertainty) of about 0.95 mA (CF), 0.3 mA (McD), 
and 0.16 mA (UHRF), which correspond to unsaturated Ca I column densities of 3.4, 1.1, and 
0.6 X 10^ cm~^, respectively. For the seven lines of sight observed both with the coude feed 
and at McDonald, the Ca I equivalent widths agree within their mutual 2a uncertainties, with no 
apparent systematic differences. The absence of absorption features at the 0.2 mA level toward a 
Vir (observed at an airmass of about 1.5) suggests that there are no significant telluric lines near 
the Ca I line. 



3. Results 

3.1. Spectra and Profile Analysis 

The normalized Ca I spectra are shown in Figure 1. For lines of sight observed both with 
the coude feed and at higher resolution, the coude feed spectrum is shown at the top, with the 
McDonald or UHRF spectrum offset below. In most cases, the vertical scale has been expanded 
by a factor of 5 or 10 to show more of the structure in these generally very weak Ca I lines. Below 
Ca I in each panel, we show the corresponding high resolution K I (or Na I Dl) and Ca II spectra, 
generally from WH, WHK, and WMH. The source of each spectrum is indicated at the right, just 
above the continuum. We have chosen the K I and Na I Dl reference wavelengths to be the /- 
value-weighted means of the wavelengths of the individual hyperfine subcomponents (7698.974 A 
and 5895.9242 A) to facilitate comparisons among the profiles. (Note the resolved Na I hyperfine 
structure, with separation ^ 1 km s~^, for several components toward e Ori; the K I hyperfine 
structure, with separation ~ 0.3 km s~^, is not seen.) 

Wc have used the method of profile fitting to determine column densities (A^), line widths (6 ~ 
FWHM/1.665), and velocities {v) for the discernible individual interstellar clouds contributing to 
the observed line profiles. We used Gaussian instrumental profiles, with the FWHM noted above, 
to fit the profiles observed with the different spectrographs. The Ca I rest wavelength (4226.728 A 
in air) and oscillator strength (1.75) were taken from the compilation of Morton (1991). In general, 
we adopted the minimum number of components needed to adequately fit the observed profiles, 
given the S /N achieved in each case — such that the rms deviation between the data and the fitted 
profile is comparable to that previously measured for the adjacent continuum. Components were 
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added, as needed, to account for apparent asymmetries and/or local systematic deviations from 
the initial simpler fits. In some cases, we used comparisons with similar high-resolution spectra of 
the stronger lines of K I, Ca II, and/or Fe I to infer more complex structure (e.g., for the main 
absorption toward (, Per and Oph) or to confirm very weak components of marginal statistical 
significance in the Ca I spectra (e.g., toward e Ori, a Cyg, and a few other stars). 

To illustrate the fitting procedure, in Figure 2 we show the individual components (dotted 
lines) and the composite fitted profiles (smooth solid lines) for two-component and four-component 
fits to the Ca I absorption observed toward o Per (solid histograms). The residuals (data minus 
fit) are shown as the lower histogram in each panel. The two-component fit to this profile (with 
velocities near 11.3 and 13.6 km s~^) is inadequate, as the rms deviation between the data and 
fit within the absorption-line region (0.0071) is noticeably larger than the value observed in the 
nearby continuum (0.0062). Moreover, the observed absorption near both those velocities appears 
asymmetric; single-component fits to the left and right halves of the profile yield rms deviations of 
0.0071 and 0.0065, respectively. Splitting either of those two components into two sub-components 
reduces the overall fit rms deviation to 0.0056 or 0.0059 — i.e., below the continuum rms deviation 
— but does not address the apparent asymmetry in the other component. Four components (at 
minimum) appear necessary, and yield a fit rms deviation of 0.0047. Comparison with the K I 
spectrum (Fig. la; observed at comparably high resolution) indicates that four components — at 
very similar velocities to those found for Ca I — are also needed to fit the strongest K I absorption 
(WH). 

The individual components derived in the fits to all the Ca I profiles are noted in Figure 1 by 
tick marks above the spectra. The derived individual component parameters are listed in Table 2, 
where successive columns give the star name, observatory, component number, heliocentric velocity, 
column density, and ^-value. Of the 112 individual interstellar cloud components listed in Table 2, 
70 were derived from UHRF and/or McDonald spectra (15 stars), with the remaining 42 from the 
lower resolution coude feed spectra (15 stars). Where we have spectra from both the coude feed 
and McDonald, the component parameters are those derived from the latter (higher resolution) 
data. 

As was found in our previous high-resolution surveys, these new higher-resolution, higher S/N 
Ca I spectra show complex component structure in most lines of sight. Where we have very high 
resolution spectra (FWHM < 0.6 km s~^) for Ca I and for the stronger lines of K I, Na I, and/or 
Ca II, the component velocity structures derived from independent fits to the various line profiles 
are generally very similar. The Ca I components generally correspond to the stronger components 
seen in K I and/or Ca II. Where we have Ca I spectra only from the coude feed (lower resolution 
and often lower S/N), comparisons with the corresponding higher resolution spectra of the stronger 
lines suggest that the true Ca I structure is more complex than we have been able to discern. 

The profile fits also yield estimates for the uncertainties in the component parameters. For 
well-defined components, the formal la uncertainties in the column densities are typically no more 
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than 10-20% (but at least several times 10^ cm^^) — comparable to the uncertainties that would 
be inferred from the error bars on the equivalent widths. The uncertainties in the column densities 
can be somewhat larger for significantly blended components, however. For most well-defined 
components in the McDonald and UHRF data, the derived velocities have formal uncertainties less 
than 0.1 km s~^, but the uncertainties are generally somewhat larger for broader or more severely 
blended components and for velocities derived from the coude feed spectra. Comparisons with 
the K I spectra suggest that the absolute velocities for most well-defined components should be 
accurate to within about 0.2-0.3 km s^^. The formal uncertainties in the line widths are typically 
of order 20% for well-defined components in the higher-resolution spectra, but can be 30-50% for 
weaker, narrower (5 < 0.5 km s~^), and/or more blended components. The nine 6-values enclosed 
in square brackets were fixed to facilitate convergence of the fits, and are about as well determined 
as the 6-values allowed to vary. The 42 6-values enclosed in parentheses, often for weak and/or 
blended components, were essentially arbitrarily (though not unreasonably) fixed in the fitting, and 
are less well determined. As for Na I, Ca II, and K I, uncertainties in the component parameters 
due to not discerning the full Ca I component structure may well be more significant than the 
formal uncertainties derived in the profile fits. 



3.2. Comparisons with Other Ca I Data 

Because the Ca I lines are typically very weak, very few high resolution (FWHM < 2 km s^^) 
spectra had previously been obtained (Hobbs 1971, using the PEPSIOS interferometric spectrom- 
eter); in many cases, only equivalent widths were reported (e.g.. White 1973; Lambert &; Banks 
1986; Gredel, van Dishoeck, &: Black 1993; Allen 1994). (Note that for the last three references, 
the primary focus was the nearby CH+ line at 4232 A, so any discussion of Ca I was rather brief.) 
In the last column of Table 1, we list equivalent widths from various prior studies for the stars 
in our current sample. In most cases, the values measured for a given line of sight are consistent 
within the listed 2a uncertainties. Some of the upper limits reported by Hobbs (1971), computed 
for single, unresolved absorption features, are lower than subsequent detections, however. In all 
such cases, at least two weak components are detected in the new spectra. 

Comparisons of the new high-resolution Ca I spectra with those shown by Hobbs (1971) and 
Fcderman &: Hobbs (1983) reveal apparent differences in several cases. (1) Toward A Ori, Federman 
& Hobbs found two components of equal equivalent width [1.9ib0.4 mA (2(t)], at heliocentric veloc- 
ities of about 6.6 and 24.2 km s^^. While the higher velocity feature is reasonably consistent with 
the main absorption shown in Figure lb, the lower velocity feature is much stronger than the weak 
component near 7.6 km s~^ seen (tentatively) in the new McDonald spectrum. There is no evidence 
for significant temporal variations in the stronger Na I D lines at that velocity, over a somewhat 
longer time period (Hobbs 1969; Watson 1998; Price et al. 2001). (2) Toward 23 Ori, roughly 30% 
of the absorption shown by Fcderman & Hobbs (of order 1.5 mA) appears to lie between about 
12 and 16 km s~^, but none is seen at those velocities in the coude feed spectrum in Figure lb. 
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While the absorption features toward A Ori and 23 Ori appear to be statistically significant in the 
spectra shown by Federman & Hobbs, we note that there are comparable, but apparently spurious 
features (both positive and negative) toward several other stars in their small sample. (3) The 
new McDonald spectrum of A Cep (Fig. le) shows a fairly strong, narrow Ca I component near 
—3 km s~^. While that component is not apparent in the spectrum shown by Hobbs (1971), that 
velocity is near the end of the PEPSIOS scan, and better agreement between the spectra might be 
obtained by adjusting the continuum in the earlier data. 



4. Discussion 

4.1. Component Statistics 

In our previous surveys of Na I, Ca II, and K I (WHK; WMH; WH), we examined the statistical 
properties of the ensembles of individual component column densities, line widths, velocities, and 
velocity separations between adjacent components (Sv). For all three species, the true median Sv 
(allowing for additional unresolved structure) is estimated to be of order 1.2 km s~^. The median 
line widths for individual Na I and K I components are both about 1.1-1.2 km (FWHM); the 
corresponding Ca II components are typically somewhat broader, however. Adjacent components 
can have rather different properties [e.g., line widths, A'^(Na l)/Af(Ca II) ratios]. Those references 
also discuss the effects of limited samples, inhomogeneities in resolution and S/N, and the likelihood 
that we have not discerned all the components actually present. 

Because the Ca I absorption is generally weak, the number of components detected in this 
survey is significantly smaller, and provides little new information on the statistics of v, Sv, and 
N; we merely note that the median individual component Ca I cohimn density is of order 2 x 
10^ cm~^. While the set of Ca I component 6- values is rather limited, however, it does provide 
interesting constraints on the temperature, internal turbulent velocity (vt), and distribution of the 
gas traced by Ca I, since b = {2kT/m + 2uf )i/^ where m is the atomic weight.^ The median h for 
the primary sample of 47 well determined 6- values derived from the high-resolution UHRF and/or 
McDonald spectra is 0.66 km s~^, which corresponds to a maximum temperature 

Txaax 1030 K or 

to a maximum turbulent velocity vt max 0.47 km s-i. Roughly 85% of the primary 6- values fall 
between 0.2 and 0.9 km s^^, corresponding to values of T^ax between 95 and 1920 K or values of 
Vt max between 0.14 and 0.64 km s^^. The median 6- values for the coude feed and combined samples 
are both somewhat larger, reflecting the less complete knowledge of the component structures in 
the lower resolution spectra. Besides thermal and turbulent broadening, unresolved blends and 
cloud stratification may also contribute to the widths of the observed Ca I lines (see WHK and 
WMH). Comparisons with the profiles of K I and Ca II do suggest that most Ca I components with 
6 > 1 km s~^ may in fact consist of unresolved blends. The "true" median Ca I 6-value therefore 
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could be somewhat smaller. 

If Ca I, Ca II (m = 40), and K I (m = 39) are similarly distributed, we would expect that 
the respective 6-values would be essentially identical, regardless of whether thermal broadening or 
turbulent broadening is dominant. For the components in the respective high resolution primary 
samples, the median 6-values are 6incd(Ca I) = 0.66 km s^^, 6mcd(Ca II) = 1.33 km s~^, and 
^med(K I) = 0.67 km s~^. We have identified 35 "corresponding" components in Ca I and K I 
for which the velocities agree to within 0.4 km s~^ and for which any immediately neighboring 
components have roughly similar properties (as for the corresponding Na I and Ca II components 
discussed by WMH). The median values of b are 0.69 and 0.65 km for the corresponding Ca I and 
K I components, respectively. The close agreement for the Ca I and K I 6-values suggests that those 
two species are essentially coextensive in the interstellar gas — and also with Na I (WH). While the 
typical uncertainties in the 6-values would allow 6(Ca I) ~ 6(Ca II) for many of the corresponding 
Ca II components with b < 1.2 km s~^, the data suggest that 6(Ca II) is often greater than 6(Ca I) 
— both globally and for individual components seen at the same velocity. Barlow et al. (1995) 
and WMH have suggested that this difference in 6-values may imply that Ca II and the three trace 
neutral species are not identically distributed, even though similarities in component velocities 
suggest that the various species are associated. For a given velocity component, Ca II may occupy 
a somewhat larger volume, characterized by a larger temperature and/or greater internal turbulent 
motions (where Ca might be less severely depleted), than do the trace neutrals. 

4.2. Comparisons of Ca I with Other Species 

In this section, we examine some of the relationships between Ca I, Ca II, K I, Na I, C I, 
S I, Fe I, H (H I + 2H2), and H2 — which (in principle) can yield insights into the abundances, 
depletions, and physical conditions in the predominantly neutral clouds in which those species 
reside. The interpretation of the observed relationships, however, involves several assumptions 
and approximations which are not strictly valid, and which can bias conclusions drawn from those 
relationships. We begin by briefly examining those assumptions and noting how they may affect 
our interpretation of the observed data. 

4.2.1. Assumptions and Approximations 

While we would like to determine abundances and physical conditions at specific locations 
within individual interstellar clouds, we do not have the data to do so. The high-resolution spectra 
available for some lines of sight suggest that most sightlines are complex, containing multiple 



*In this definition of 6, vt is the one-dimensional rms value for a Gaussian spectrum of internal turbulent velocities; 
the full three-dimensional turbulent velocity is thus 
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blended components (clouds) which can be characterized by different relative abundances and 
physical properties. For example, it is clear from the profiles in Figure 1 that Ca II is often detected 
(and sometimes strong) over a much wider velocity range than Ca I and K I. For many species and 
many sightlines, however, we have only total column densities integrated over the whole line of sight 
— whether because of low spectral resolution (for most UV lines) or, for H I and H2, the inextricable 
blending of very strong, broad individual components. Where we do have very high-resolution 
spectra, we may obtain more detailed information by comparing column densities for corresponding 
individual components identified in fits to the line profiles (e.g., Welty et al. 1999b) or, alternatively, 
by comparing "apparent column densities" as functions of velocity (e.g., Jenkins &: Tripp 2001). 
Even then, however, it is likely that individual clouds are not homogeneous, and that different 
species may be distributed somewhat differently within the clouds. For example, considerations of 
ionization equilibrium suggest that trace neutral species will be more concentrated in the denser 
parts of the clouds than the corresponding dominant first ions; differences in photoionization cross- 
section and depletion behavior may lead to further differences in distribution among the various 
trace and dominant species. The generally larger 6-valucs found for Ca II, compared to those found 
for Ca I and K I, may reflect a combination of those effects. 

It is commonly assumed that photoionization and radiative recombination dominate the ion- 
ization equilibrium of heavy elements in diffuse clouds, such that 

r(XO)n(XO) = a,(XO,r) nen(X+), (1) 

where F is the photoionization rate and is the radiative recombination rate coefficient. It is also 
usual to assume (despite the potential problems noted above) that the local volume densities (n) 
may be replaced by the corresponding column densities (N) — either for individual components 

or (in most cases) for the total line of sight. Ratios of trace neutral species with the corresponding 
dominant first ions [measured directly or estimated from A^(H) and assumed depletions] thus can 
yield estimates for rie (if the temperature and radiation field can be constrained) and the local 
hydrogen density nn (if the fractional ionization ne/nn is assumed to be equal to the gas phase 
abundance of carbon). 

In the discussions below, we will use the equations of photoionization equilibrium, with the 
photoionization rates calculated for the WJl radiation field (de Boer et al. 1973) and the radiative 
recombination rates at T = 100 K (usually), as tabulated by Pequignot & Aldrovandi (1986). We 

define the electron density derived under the assumption that radiative recombination dominates 
as ne(rad) = (F/a,.)[A^(X I)/A^(X II)]. Unless otherwise noted, we use total line-of-sight column 
densities, since individual component values generally are not available for hydrogen, C I, or S I. 
Column densities for Ca I, Ca II, Fe I, and S I (with references and a brief discussion of uncertainties) 
are given in Appendix B of this paper (Table 9); values for the other species are given in the appendix 
to WH. Wc also adopt the solar system meteoritic reference abundances (Aq) listed by Anders & 
Grevesse (1989), except for carbon, for which we take the value for the solar photosphere from 
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Grevesse k, Noels (1993)^. The reference abundances and rates are collected in Table 3. 

While we will adopt those assumptions (at least initially) in the discussions that follow, we will 
find in some cases that they lead to results that are inconsistent with those implied by other diagnos- 
tics of the physical conditions within the clouds — suggesting that some (or all) of the assumptions 
may not be valid. For example, the electron densities inferred from different A^(X l)/A^(X II) ratios 
in a given sightline can differ by an order of magnitude or more (e.g., Fitzpatrick & Spitzer 1997; 
Welty et al. 1999b). Moreover, the estimates for Ue and tzr can be systematically higher than the 
values derived from analyses of the fine-structure excitation of C II and C I (Welty et al. 1999b; 
WH). Such systematic differences suggest that other processes besides photoionization and radia- 
tive recombination may have a significant effect on the ionization balance of heavy elements in at 
least some diffuse clouds. 

One promising candidate for enhancing the various trace neutral species is charge exchange 
between the corresponding singly ionized species and neutral and/or negatively charged large 
molecules/small grains (Lepp et al. 1988; WH; Weingartner & Draine 2001; Liszt 2003). If this 
"grain-assisted" recombination dominates over radiative recombination, we have 

r(XO) n(XO) = a5(X0,V,r) nn n(X+), (2) 

where ag is the recombination rate coefficient {per hydrogen nucleus). In this definition, ag includes 
the effects of the grain abundance, the grain size distribution, and the grain charge distribution; Ug 

depends on rig, T, and the radiation field through the "charging parameter" ^ = G\/T Inf.., where 
the energy density in the ambient interstellar radiation field between 6 and 13.6 eV is G times that 
of the Habing (1968) field. In the next-to-last column of Table 3, we list values of ag for the "typical 
cold, neutral medium" (T = 100 K, ^ = 400 K^/^ cm^), calculated via equation 8 in Weingartner 
& Draine (2001). 

Equation 16 of Weingartner k, Draine (2001) gives the electron density for the more general 
case in which both radiative and grain-assisted recombination may be significant: 



R r 
ne = 

Or 



^ _ /I - s + i?\ CKp nn 

~ V 1 + ii: 7 R r 



(3) 



where R = N(X l)/N(X II) is usually much less than 1.0, s is the probability that an ion sticks 
to the grain after charge exchange, and the quantity in front of the square braces is ne(rad). For 
any assumed s, the effect of including grain-assisted recombination thus depends in a somewhat 
complicated way on the column density ratio R, the local density nn, the temperature (mildly), 
and the strength of the radiation field — directly and/or through the dependence of Og on those 
factors. In general, the effect is largest for small R (and thus usually for small inferred rig), for large 
nn, and for small T/ng (e.g., for Ca^, Na, Mg, and K). For s = 1 (i.e., all the neutralized ions are 



^Recent slight revisions for some of the reference abundances (e.g., Grevesse & Sauval 1998; Holweger 2001) would 
have a minor impact on the results discussed in this paper. 
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incorporated into the grains, and not returned to the gas phase), the quantity in the square braces 
becomes [1 — ag nn / F] ~ 1 (unless tt-h is much greater than 10^ cm~^), so that ng ~ rig (rad). For 
s = (which maximizes the effect of the grain-assisted recombination), the above equation can be 
re-arranged to yield R = ne(rad) [1 + (ug nH)/(ar n-g)]. The ratio {ag nu)/{ar Wg), effectively a 
measure of the relative contributions of grain-assisted and radiative recombination, can range from 
about 0.1-10 for fractional ionizations between 2 x 10~^ and 0.01, for the recombination coefficients 
listed in Table 3. The effect of grain-assisted recombination thus can vary for different elements in 
any given line of sight (depending on ag/ar and s), and can also vary from one sightline to another 
for any given element (depending on R, T, F, and nn)- 

In the following sections, we examine a number of possible relationships between individual 
ions and also between ratios of different atomic and molecular species, in an attempt to determine 
which physical processes have an important effect on the abundances of the various species. In 
each case, we first determine whether or not a correlation exists between the two quantities, using 
several different estimates for the correlation probability. The first estimate is that associated with 
the linear correlation coefficient, which is derived from the values of all points detected in both 
variables (e.g., Bcvington 1969). As we have chosen to weight the data (using l/o"^), the tabulated 
linear correlation coefficients and corresponding probabilities arc averages of the values determined 
for considering each variable as the dependent one. Because only limits (for one or both of the 
variables) are available for many of the points in some of the possible correlations, we also provide 
estimates for the correlation probabilities using the three bivariate correlation tests in the survival 
analysis (or censored data analysis) package ASURV (Isobe, Feigelson, & Nelson 1986), which use 
both detections and limits.^ (Note that in all cases we give the probabilities of no correlation 
— i.e., a small value for the probability corresponds to a strong correlation.) While the four 
different methods yield concordant results for most of the possible correlations, there are several 
cases for which the survival analysis tests suggest that a correlation is present, even though the 
linear correlation coefficient suggested otherwise. For the confirmed correlations, we then perform 
a regression analysis which allows for uncertainties in both variables (but which uses only the 
detected points), and which yields both the slope and intercept of the "typical" relationship and 
an estimate of the scatter (WH).^ For some of the relationships, we also perform the fits with the 
slope fixed at 1.0 or 0.0 for easier comparison with theoretical approximations. 



® Specifically, we used Cox's proportional hazard model, the generalized Kendall's tau method, and the generalized 

Spearman's rank order method. Version 1.2 of the ASURV package (LaValley, Isobe, & Feigelson 1992) was obtained 
from the statistical software archive maintained at Penn State (http://www.astro.psu.edu/statcodes). These routines 
do not allow weighting of the data. 

'^We used a slightly modified version of the subroutine regrwt.f, obtained from the Penn State statistical software 
archive, to perform the regression fits. While WH assumed equal weights, here we weight the data by l/cr'^; the 
weighted and unweighted fits generally yield similar results, however. The bivariate regression routines in survival 
analysis packages like ASURV do not use uncertainties (for either variable), and cannot be used when there are both 
upper and lower limits for either variable. 
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As much as possible, we attempt to minimize biases resulting from differences in the spatial 
distribution of different species (as noted above) by comparing or considering ratios of species that 
are likely to be similarly distributed, by virtue of having similar dependences on temperature, 
density, and/or ambient radiation field. For example, both considerations of ionization equilibrium 
and empirical evidence suggest that the various trace neutral species (e.g., C I, Na I, K I) are 
fairly well correlated with each other — though there may be some differences in ionization and/or 
depiction (as discussed below). Moreover, the column densities of both Na I and K I exhibit fairly 
well-defined relationships with the column densities of the molecular species CH and H2 (with 
some exceptions; WH). Most dominant species (e.g., Si II, Fe II, Zn II) should also be fairly well 
correlated with each other, though (again) there will be differences due to depletion behavior (see, 
for example, the absorption-line profiles and curves of growth for the line of sight to 23 Ori in Welty 
et al. 1999b). Even Ca II, whose behavior exhibits an interplay between ionization and depletion 
effects, shows a surprisingly good correlation with Fe II, which is typically a "dominant-depleted" 
species in the primarily neutral clouds sampled in this study (Welty et al. 1996, 1999a, 1999b). 

For example, where column densities for several trace neutral species are available, we may 
estimate relative abundances/depletions for those elements by taking ratios of the ionization equi- 
librium equation: 

'^(X)Ae(X) jV(X II) ^ (r/a)x NjX I) 

,5(Y)Ao(Y) ~ N{Y II) (r/a)Y A^(Y I) ' ^ ^ 

where S{X) is the depletion of X, Aq is the solar reference abundance of X, and we have assumed 
that the singly ionized species arc dominant. In taking such ratios, the dependences on the density 
(rig and/or njj), on the temperature [through a(T)] and on the radiation field (through F) — 
all of which can be rather uncertain arc substantially reduced (York 1980; Snow 1984). If 
radiative recombination dominates the production of the trace neutral species, then the effective 
recombination rate coefficient a would be equal to a^; if charge exchange with grains dominates, 
then a ~ CKg. If both processes are significant, then a would be replaced by ar[l + {ag nH)/{ar 
Ue)] or (equivalently) by ag[l + (a^ ne)/{ctg ^^H)] (for s = 0; see above). The ratios n(X Il)/n(Y II) 
and n(X l)/n(Y I) should vary much less through each cloud than either the individual densities or 
the trace/dominant ratios. Because the depletion of carbon in the local ISM is apparently uniform 
[(5(C) ~ 0.4, relative to solar abundances (Cardelli et al. 1996; Sofia et al. 1997)], the roughly 
constant ratios of the column densities of Li I, Na I, and K I, relative to that of C I, have thus 
yielded values for the typical depletions of lithium, sodium, and potassium: D(X) = log (5(X) ~ 
—0.6, —0.6, and —0.7 dex, respectively, for a ~ (WH). 



4-2.2. Calcium Ionization Balance 

Measurements for two adjacent ionization states (Ca I and Ca II) yield constraints on the 
ionization of calcium in interstellar clouds. Because the ionization potential of Ca II is about 11.87 
eV, Ca III is expected to be the dominant form of calcium in very diff'use clouds. The combination of 
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a relatively small photoionization rate for Ca II and a relatively large radiative recombination rate 
coefficient for Ca III (Table 3), however, suggests that Ca II could be dominant in somewhat thicker 
clouds. Taking the equilibrium equations for two successive ionization states, A^(Ca l)/A^(Ca II) 
= ne/(r/Q^)cai and iV(Ca Il)/iV(Ca III) = ne/(r/a^)caii, we obtain N{Ca Il)/iV(Ca III) ~ 1300 
A^(Ca I) /iV(Ca II). For individual components with A^(Ca I) /iV(Ca II) in the observed range 0.0015- 
0.08, we thus have iV(Ca Ii)/N{Ca III) ~ 2-100 — so that Ca II would be the dominant form of 
calcium wherever Ca I is detected. 

Several additional factors may modify that result, however. (1) If Ca II is more broadly 
distributed than Ca I within individual clouds (as discussed above), then the Ar(Ca II)/Ar(Ca III) 
ratio would be even larger in the region occupied by Ca I. The magnitude of this effect is difficult to 
estimate, however. (2) On the other hand, if the recombination is dominated by charge exchange 
with small grains [or large molecules, with ag{Ca III) ~ 3 x ag{Ca II) (Weingartner &: Draine 
2001)], then A^(Ca II)/iV(Ca III) would be reduced by a factor of about 1.7, and Ca III would be 
comparable to Ca II in cases where N{Ca I)/N{Ca II) is low. (3) Moreover, if the unattenuated 
{T/ar)cai is actually smaller than the value listed in Table 3 (a possibility noted below), then 
A^(Ca Il)/A^(Ca III) would be correspondingly reduced, implying more (unobserved) Ca III, and 
thus less severe depletion of calcium than inferred from Ca II alone [again, where A''(Ca l)/iV(Ca II) 
is low]. 



4-2.3. Ca I (and Other Trace Neutral Species) vs. K I 

Welty k. Hobbs (2001) found fairly tight, essentially linear relationships between the logarith- 
mic column densities of Li I, Na I, and K I, with slopes near unity and rms scatter less than about 
0.2 dex. The close correspondence between those three species suggests that they respond very 
similarly to changes in physical conditions in the ISM and that they are generally coextensive — 
which is not surprising, in view of similarities in ionization cross-section and chemical behavior. 
(Slight differences were seen, particularly for Na I, for sightlines in the Sco-Oph region and in 
several other regions, however.) While the column density of C I also showed a strong correlation 
with the column densities of Na I and K I, the slightly steeper slopes (1.1-1.3) may reflect the 
differential effects of extinction on the various photoionization rates in the higher column density 
sightlines (see below) and/or uncertainties in the C I /-values (Jenkins k, Shaya 1979; WH; Jenkins 
& Tripp 2001). 

The additional data collected in Appendix B of this paper permit similar comparisons between 
A^(K I) and the column densities of Ca I, S I, and Fc I. For example, in Figure 3 we plot the 
total sightline values of log[A^(Ca I)] vs. log[A^(K I)]. While the correlation probabilities listed in 
Table 4 indicate strong correlations in all three cases, the fitted regression slopes (Table 5) all differ 
from unity. The slope of the best-fit (solid) line for log[iV(Ca I)] vs. log[iV(K I)] is 0.60±0.06 
(Table 5; all sightlines with detections of both species). Depending on the sample used (including 
the Sco-Oph sightlines or not), the slopes for log[A?^(S I)] and log[A?^(Fe I)] are about 1.2-1.5 and 
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0.8-0.9, respectively. The rms scatter for the unconstrained fits, in all cases between 0.12 and 0.21 
dex, is somewhat larger than the values expected from the nominal uncertainties on the individual 
column densities. Either the uncertainties have been underestimated or there are real deviations 
from the mean relationships in some sightlines. The slight reductions in the scatter when the Sco- 
Oph sightlines are excluded and the somewhat larger scatter for the correlations involving Ca I 
and/or Fe I (which are expected to be more affected by variations in depletion) both suggest that 
real deviations are present. 

In principle, several effects could lead to slopes greater than or less than unity in comparisons 
between log[A''(X I)] and log[iV(K I)]. One possibility is that of differential extinction effects (noted 
above for C I). Differences in photoionization cross-section for the various neutral species, coupled 
with the wavelength dependences of both the interstellar radiation field and the extinction due to 
dust, imply that ratios of the photoionization rates [r(X'^)/r(Y'^)] will vary somewhat with depth 
within individual interstellar clouds (e.g., Roberge, Dalgarno, &; Flannery 1981). In column 5 of 
Table 3, we give the change in log(r) between the edge (the unattenuated value) and the center of 
a cloud with total visual extinction 1 mag, for grain model 2 of Roberge et al. (drawn from their 
Figure 2b). Those values suggest that C I, S I, and Ca II (with the highest ionization potentials) 
would be slightly enhanced (by 0.15-0.24 dex), relative to K I, in going from the edge to the center 
of the cloud, while species such as Mg I, Al I, and Si I would be slightly reduced (by 0.10-0.14 
dex). The relative reductions of Na I, Ca I, and Fe I would be even smaller (< 0.1 dex). The 
relationships between N{K I), A?^(H), and E{B — V) discussed by WH suggest, however, that even 
the thickest (highest column density) individual clouds in our sample probably have total visual 
extinctions of only about 0.6 mag — so that the effects of differential extinction would be even 
smaller than those estimated above (especially when integrated over the whole cloud). While the 
slightly steeper slopes found for C I and S I thus might be due (at least in part) to differential 
extinction in the thicker clouds, it seems unlikely that the shallower slopes found for Ca I and Fe I 
could be explained in that way. 

A second possibility is that of differences in how the recombination rates of the various singly 
ionized species might vary — either with depth within an individual cloud or between more diffuse 
and denser clouds. As an example, we consider the ratio A^(S I)/A^(K I), including radiative and/or 
grain-assisted recombination, for T = 100 K and tp = 400 K^/^ cm"^ (the "typical" values) and for 
T = 25 K and = 100 K^^"^ cm'^ (representing colder, denser regions). We choose those two species 
for several reasons: (1) differences in depletion are unlikely to be significant, (2) there are (slight) 
differences in the temperature dependence of (Pequignot & Aldrovandi 1986), (3) there are 
differences in the ratio ar/ag, and (4) the observed relationship is steeper than linear. If radiative 
recombination dominates, the ratio ar{S)/ar(K) is equal to 1.91 for the "typical" case and to 1.69 
for the "colder, denser" difference of 13% (0.05 dex). If grain-assisted recombination 

dominates, the ratio ag{S)/ag(K) equals 1.36 and 1.44 for the two cases — a difference of 6% (0.02 
dex). When both radiative and grain-assisted recombination arc significant, we take ratios of the 
quantity ag[l + (a^ ne)/{ag nn)] (for s = 0; see § 4.2.1 and § 4.3.1) for S I with respect to its value 
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for K I; we further assume a constant fractional ionization ng/nn = 5 x 10~^. For the "typical" 
case, the ratio is equal to 1.50; for the "colder, denser" case, the ratio is equal to 1.53 — only 
2% larger. If the fractional ionization decreases at higher densities (a possibility noted below), the 
"colder, denser" ratio will decrease slightly, to a minimum of 1.44 (the ratio of the ag). These 
rather slight differences suggest that differences in recombination rate behavior are not responsible 
for the observed steeper than linear relationship between N{S I) and N{K I). 

A third possibility is that of differences in depletion behavior (e.g., with depth or density) 
for the different elements. Observed increases in the severity of depletions (e.g., for iron, silicon, 
magnesium, calcium) for sightlincs with higher mean densities (< nn >) or higher fractions of 
hydrogen in molecular form /(H2) have suggested that the depletions may depend on the local 
hydrogen density nn (e.g., Jenkins 1987; Cardelli 1994). Wakkcr & Mathis (2000) have shown 
similar trends for the depletions of several elements versus A^(H). Given the roughly quadratic 
dependence of A'^(K I) on A^(H) and the apparent lack of significant variation in the potassium 
depletion (WH), we thus would expect (1) that the ratio N{X 1)/N{K I) — which is roughly 
proportional to the depletion of X (sec next section) — would decrease at higher overall column 
densities, and (2) that the slope of the relationship between log[A^(X I)] and log[A^(K I)] would 
be less than 1.0 for elements X which can be severely depleted. For example, Wakker & Mathis 
find that N{Fe Il)/iV(H) is roughly proportional to [Af(H)]~0-6, for lO^*^-^ cm'^ < iV(H) < lO^^-^ 
cm-2. If N(K I) is proportional to [iV(H)]2, then N{Fe l)/iV(K I) oc N(Fe Il)/iV(K II) would 
be roughly proportional to [N{K I)]~°'^, so that N{Fe I) would be proportional to [N{K 1)]°"^ — 
slightly weaker than the observed dependence. Differences in depletion behavior thus can account 
quite naturally for the shallower slopes found for N[Ca I) and A^(Fe I) [and also for Mg I, for which 
fewer data are available], versus A''(K 1). 



4.24. Cal/K I and Ca I/Ca II vs. f (H2) 

As noted above, the components detected in Ca I generally have counterparts in both K I 
and Ca 11. For Ca 1 and K I, the generally good correspondence in both velocity extent (Figure 1) 
and component 6- values (§ 4.1) suggests that the two species are essentially coextensive. The 
typically larger 6- values for Ca II, however, suggest that Ca II is somewhat more broadly distributed; 
moreover, Ca II is often detected over a wider velocity range. For the set of individual components 
detected in all three species, the column density ratios A^(Ca l)/A^(K I) and A^(Ca l)/A^(Ca II) 
can range over factors of about 500 and 50, respectively — though the ranges for components 
within any individual sightline are generally much smaller than those extreme values. The total 
line-of-sight values (or limits) for A/'(Ca l)/N{Ca II) will underestimate the "true" values (i.e., for 
the components detected in Ca I) for cases where Ca II is detected over a much wider velocity 
range — but generally not by more than a factor of 2. In this section, we use the equations of 
photoionization equilibrium to interpret the behavior of those two column density ratios in the 
different environments probed by the sightlines in our sample. 
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In photoionization equilibrium, the ratio 

A^(Ca I) {V/ar)K [(5(Ca)A0(Ca)] 
iV(KI) ~ (r/a,)ca [5{K)Kq{K) 



12.4 5(Ca) (5) 



is sensitive primarily to the calcium depletion, with only weak dependence on T or the radiation 
field. We have assumed a constant potassium depletion [5(K) = 0.2] and that Ca II and K II are 

the dominant ionization states. The ratio A'^(Ca l)/A^(Ca II) = ne/{T/ar) depends on the electron 
density, the ambient radiation field, and the temperature (since ~ cto T~^, with /3 typically 
0.6-0.7). If the fractional ionization rie/n^i is assumed constant, then the latter ratio would also 
depend on nn- If charge exchange with grains were to dominate over radiative recombination, 
then N{Ca l)/N{K I) ~ 7.6 S{Ca) [with 5(K) = 0.36] and iV(Ca I)/iV(Ca II) = nu/ir/ag). If 
the sticking parameter s ~ 1 for calcium (typically severely depleted), however, then radiative 
recombination would be responsible for most of the Ca I, and we could have 

iV(Cal) (r/a,)K He [(5(Ca)Ao(Ca)] 3 

W(Kiy~(r7^^ TOAo(K) 

In Figure 4, we compare the two ratios A^(Ca l)/A^(K I) and A^(Ca l)/A^(Ca II) with the 
fraction of hydrogen in molecular form /(H2), which is often considered to be indicative of the 
local density nn (e.g., Cardelli 1994). In this figure, the filled circles represent the total line-of- 
sight values, while the crosses represent individual components (only those detected in the relevant 
species). Sightlines in the Sco-Oph and Orion regions, some of which can exhibit deviations from 
the otherwise well-defined mean relationships in some comparisons (WH), arc marked with open 
circles and open squares, respectively. The individual components are assumed to have the same 
/(H2) as the sightline as a whole — which is likely a reasonable approximation for sightlines with 
low overall /(H2), but not for minor components in sightlines with high /(H2). For higher densities, 
we might expect lower temperatures, higher molecular fractions, and higher ne/{T/ar) ratios. For 
stronger radiation fields, we expect lower /(H2). 

The top plot, log[A^(Ca I)/A^(K I)] versus log[/(H2)], indicates that A^(Ca I)/A^(K I) is roughly 
constant for /(II2) < 0.1, but decreases markedly for /(H2) > 0.1. This behavior is very similar 
to that found for the gas phase abundances of dominant ions of elements whose depletions become 
more severe in colder, denser clouds [e.g., A^(Fe II)/A^(H)] — consistent with the expectation that 
N{Ca I)/A^(K I) is an indicator of the depletion of calcium. 

In the bottom plot, however, there is no apparent dependence of A^(Ca l)/A^(Ca II) on /(H2) 
— suggesting cither that ne/{T/ar) does not depend on the local density (see also § 4.3.3 below) 
or that /(II2) is not a very reliable indicator of nn or that the assumption of photoionization 
equilibrium does not hold. The apparent lack of dependence on /(H2) is particularly interesting 
for sightlines with relatively high molecular fraction (log[/(H2)] > —1.5) — where we expect that 
the H2 lines will be strong and self-shielded and where we do see variations in depletion. 
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4-2.5. Is Ca I Enhanced Relative to Other Trace Neutrals? 

In the main clouds toward 23 Ori, the electron density inferred from the ratio A'^(Ca I) /iV(Ca II) 
is much higher than the values derived from other such ratios (Welty et al. 1999b; see below) — 
perhaps suggesting that some as yet undetermined process acts to enhance Ca I, relative to the 
other trace neutral species. In this section, we examine this issue for a larger sample of sightlines. 

In photoionization equilibrium, we have 

iV(Cal) (r/a^)x 1 A^(Call) 

N{X I) ~ (r/a,)ca S{X)Aq{X) 7V(H) 

For the elements X = C, Na, K, and S, we would expect essentially linear relationships between 
A'^(Ca l)/N{X I) and A^(Ca Il)/iV(H), with only weak dependences on the temperature or radiation 
field, as the depletions of those four elements are thought to be uniform (within factors of 2). In 
Figure 5, we plot Iog[iV(Ca 1)/N{X I)] versus log[A'"(Ca II)/A^(H)] for those four elements. The 
solid lines of unit slope correspond to the relationships expected for solar (undepleted) abundances; 
the dotted lines correspond to the relationships expected for typical depletions of the four elements. 
If grain-assisted recombination were much more important than radiative recombination for both 
calcium and the other four elements, the solid and dotted lines would shift vertically by less than 
about 0.2 dex in all cases. If radiative recombination dominates for calcium but grain-assisted 
recombination dominates for the other four elements, then we would have 

iV(Cal) {T/ag)x 1 iV(Ca II) 

7V(X I) ~ (r/a,)ca ^ 5(X)A0(X) iV(H) 

While the correlation probabilities collected in Table 4 suggest that correlations are present 
in all four cases, the slopes of the relationships are rather uncertain. The regression fits to the 
detected points (weighted; denoted by an R in the last column of Table 5) yield slopes ranging 
from 1.45 to 2.25, but the fits which attempt to account also for the limits (unweighted; denoted 
by an A)^ yield slopes from 0.5 to 0.9. For comparison with the "expected" linear relationships, 
the dashed lines in Figure 5 show the fits to the detected points, with the slopes fixed at 1.0. The 
observed data generally lie above the expected (depleted) lines, with mean offsets of factors of 
3.2-6.9 — as if A''(Ca I) were systematically higher than expected. Because calcium is generally 
more severely depleted in the higher Ar(H) components (where the trace neutral species are more 
abundant), the 7V(Ca II)/A?^(H) ratio is likely to be lower there than for the line of sight as a whole 
— making the discrepancies between the actual and predicted A^(Ca I)/A^(X I) ratios even larger. 
In all four cases, the Sco-Oph sightlines stand apart, with A^(Ca l)/A^(X I) higher than for other 
sightlines at comparable Af(Ca Il)/Af(H) [or Af(Ca Il)/iV(H) lower for a given iV(Ca I)/A^(X I)]. 
That apparent regional difference and the overall considerable scatter — ranges in A'^(Ca I)/A'^(X I) 



*We used the EM regression algorithm in ASURV, after removing any lower limits; the Buckley-James method 
yielded similar results. 



(7) 



(8) 
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of an order of magnitude or so at any given N{Ca II)/Ar(H) — both suggest that there are real 
variations in the ratios that may reflect difl'erences in environment (e.g., through differences in the 
relative contributions of radiative and grain-assisted recombination) . 

It is difficult to understand why Ca I should be enhanced, relative to other trace neutral 
species, in many diffuse clouds. As noted above, the generally small Ca I 6-values imply that any 
selective enhancement of Ca I cannot be due to dielectronic recombination in warm (T > 3000 
K) gas. Dissociative recombination of CaH+ does not appear to be likely, either (Weingartner & 
Draine 2001). If the ratio T/a^ has been overestimated for Ca^ (by perhaps a factor of 4-5), then 
the predicted lines would lie below the data points in Figure 5, as observed. A spuriously high 
T/ar would also yield systematically high electron densities from iV(Ca l)/N{Ca II), as found for a 
number of sightlines (see next section). As already noted above (§ 4.2.2), a much smaller true T/ar 
would have potentially significant implications for both the ionization balance and the depletion of 
calcium where the N{Ca l)/N{Ca. II) ratio is small. There are, however, some sightlines where Ca I 
does not appear to be enhanced (e.g., ( Oph; see below). Some of the scatter (larger than expected 
from the nominal uncertainties) and possible regional variations seen in Figure 5 might be a result 
of differences in far-UV extinction and/or radiation fleld (e.g., for some of the Sco-Oph sightlines), 
given the difl'erences in ionization potential and cross-section among the various neutral species 
(e.g., Jenkins & Shaya 1979; Roberge et al. 1981; Welty 1989; WH). It seems likely, however, that 
several (as yet undetermined) processes may contribute both to the systematic offsets and to the 
range of observed variations seen in Figure 5. 



4.3. Electron Density and Fractional Ionization 

4-3.1. Relative Ue from Different X I/X II 

Several recent studies have suggested that the column density ratios A^(X l)/A^(X II), for 
different elements X but for the same assumed temperature and radiation field, can yield markedly 
different values for ng. For various components toward the halo star HD 215733, Fitzpatrick & 
Spitzer (1997) found systematic offsets between the electron densities derived from carbon (-1-0.4 
dex) and magnesium (+0.15 dex; using more recent /-values for the Mg II A1239 doublet), relative 
to the values obtained from sulfur. For the highest column density clouds toward 23 Ori, Welty et 
al. (1999b) found that the average rif, derived independently from nine pairs of neutral and singly 
ionized species ranged from 0.04-0.05 cm"''^ (S, Mn, Fe) to 0.09-0.26 cm"^ (C, Na, Mg, Si, K) 
to 0.95 cm~^ (Ca). After considering various possible explanations for those differences, Welty et 
al. conjectured that charge exchange between protons and the trace neutral species might account 
for the lower rig derived from sulfur, manganese, and iron, based on rather crude estimates for 
the relevant rates given by Pcquignot &: Aldrovandi (1986). No plausible explanation was found, 
however, for the high inferred from iV(Ca l)/A'^(Ca II). Similar differences have been found for 
a few elements in the main absorption components toward the more heavily reddened stars HD 
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192639 and HD 185418 (Sonnentrucker et al. 2002 and in prep.). 

In view of those observational results, Weingartner &; Draine (2001) have considered whether 
charge exchange between singly ionized atoms and neutral and/or negatively charged small grains 
(or large molecules such as polycyclic aromatic hydrocarbons) might play a significant role in 
the ionization balance of various elements (e.g., Lepp et al. 1988). For plausible abundances of 
the small grains and models for the grain charging, Weingartner &: Draine find that this "grain- 
assisted recombination" can be comparable to radiative recombination in cold, diffuse clouds. The 
corresponding recombination coefficients (a^) per hydrogen atom for such clouds, listed in the next- 
to-last column of Table 3, range from 0.75-2.34 x 10~^^ cm^ s~^. While that small range in ag 
would seem to suggest that this process is no more element-specific than radiative recombination 
(compare the values of in column 6 of Table 3), significant differences can arise if different 
elements are characterized by very different sticking parameters (s). This additional source of 
neutral species does reduce the electron densities inferred from some of the A'"(X I)/A^(X II) ratios, 
but it does not seem able to reconcile all the element-to-element differences in Ue- In order to 
obtain more consistent rig, Weingartner &; Draine argue that dissociative recombination of CH"*" 
may enhance C I and that sodium and potassium could be essentially undepleted.^ With the 
additional assumption that s = for all elements, they thus find Ue ~ 0.03 cm^^ for carbon, 
sodium, potassium, manganese, and iron in the main clouds toward 23 Ori, consistent with limits 
derived from the fine-structure excitation equilibrium of C II and Si II. Significant, unexplained 
differences remain, however, for the (higher) electron densities derived from magnesium, silicon, 
and (especially) calcium. 

Because a number of the trace neutral species are significantly stronger than usual (relative to 
hydrogen) toward 23 Ori, one might ask whether the trends in ng found in that line of sight are 
representative of interstellar clouds in general. In Table 6, we list column densities for a number of 
neutral and singly ionized species, based on high-resolution optical and UV spectra, for the main 
interstellar components toward 23 Ori and four other stars with fairly extensive data for the trace 
neutral species. The UV spectra — generally obtained with the iifST/GHRS — were retrieved 
from the STScI archive and were fitted using profile structures determined from higher resolution 
optical spectra (as for 23 Ori). Values in parentheses for some of the dominant first ions have been 
estimated from N{Zii II) and typical depletions (Welty et al. 1999b; WH). The electron densities 
inferred from the ratios A'^(X l)/iV(X II) are given on the third line in each set of values. In each 



®Welty et al. (1999b) had assumed depletions of —0.5 dex for both sodium and potassium — slightly less severe 
than the values derived by WH from comparisons of C I, Na I, and K I, using the observed (apparently constant) 
depletion of carbon and the assumption of photoionization equilibrium. If grain-assisted recombination (with rates 
ag) were to dominate over radiative recombination, then those same comparisons would imply that both sodium and 
potassium typically are depleted by about —0.45 dex. If, however, the more stringent criterion for electron transfer 
sketched by Weingartner & Draine is adopted, then the depletions of both sodium and potassium could be even less 
severe, as the modified rates (ag) for those two elements appear to be more significantly reduced in cold, neutral 
clouds (see their Fig. 2). In principle, other possible processes, such as dissociative recombination of 011+, would 
also affect these depletion estimates. 
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case, the first entry was derived assuming photoionization equilibrium at T = 100 K, with the 
standard WJl radiation field (multiplied by 5.0 for 1 Sco and 5 Sco; see § 4.3.4). The left-hand 
panel of Figure 6 shows those rig [= ne(rad)] for the five lines of sight, relative to the respective 
averages (ng) of the values for carbon, sodium, and potassium [which are not entirely independent, 
as the typical depletions of sodium and potassium were estimated via the observed correlations 
between the corresponding trace neutral species (WH)]. In effect, the relative electron densities are 
given by 

He {r/ar)x N{X I) iV(CNaK II) 

no ^ (r/a,)cNaK iV(CNaK I) iV(X II) ' ^ ^ 

where the values for "CNaK" represent suitable averages for carbon, sodium, and potassium; we 
have assumed (for the present) that radiative recombination dominates. These relative ng should 
be fairly insensitive to the actual values of T and the ambient radiation fields characterizing the 
clouds in the different lines of sight, as well as to differences in distribution between the neutral 
and singly ionized species (§ 4.2.1; § 4.2.3). 

In order to gauge the effects of grain-assisted recombination on the inferred electron densities 
(for each clement in each line of sight), we solved equation 16 of Wcingartner & Draine (2001) 
(given above, in § 4.2.1). We have assumed that the depletions of sodium and potassium are less 
severe by about 0.2 dex (as would be the case if grain-assisted recombination dominates); values 
of nn are given for these five lines of sight in Table 7. The second Ue entries for each element in 
Table 6 give the values obtained assuming that both radiative and grain-assisted recombination 
must be considered; the right-hand panel of Figure 6 shows the corresponding relative rig (again 
normalized to the average of the values for carbon, sodium, and potassium). The entries in the 
table and the circles in the figure correspond to our assumed sticking parameters: s = for carbon, 
sodium, sulfur, and potassium (little depleted); s = 0.5 for magnesium, silicon, and manganese 
(moderately depleted); and s = 1 for calcium and iron (heavily depleted). (Weingartner & Draine 
did not consider lithium, aluminum, chromium, or nickel.) The associated "error bars" (attached 
via vertical dotted lines) show the range in relative rig if we assume either s = or s = 1 for all 
elements. In principle, those two extremes give the maximum and minimum effect of grain-assisted 
recombination; the circles in Figure 6 do not always lie between the values for s = and 1 because 
we have assumed different values of s for different elements and because we have normalized by the 
mean Ue for carbon, sodium, and potassium (which all have s = 0) . The most significant reductions 
in Ue generally are seen for carbon, sodium, sulfur, and potassium (all with ,s = 0); the less severe 
depletions for sodium and potassium also play a role. For example, the electron densities inferred 
from those four elements for the main clouds toward C Oph are smaller by factors of 3-12 than the 
values obtained by assTiniing only radiative recombination; somewhat smaller effects are seen for 
the other four lines of sight. For CaP, the effect of grain-assisted recombination is generally very 
minor, as the "radiative" ng are typically fairly large and (especially) as we have adopted s = 1. 

It is apparent from Figure 6 that the differences in relative rig found toward 23 Ori are neither 
unique nor universal. The "radiative" electron densities inferred for the main clouds toward ( Ori, 
for example, are significantly lower than those toward 23 Ori, but show a very similar pattern — 
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with roughly similar values from carbon, sodium, magnesium, and potassium (0.03-0.07 cm~^), low 
values from sulfur and iron (<0.008 cm~^), and a high value from calcium (0.30 cm~^). The patterns 
seen toward 1 Sco and 8 Sco are also similar, except that the electron densities derived from silicon 
and iron are comparable to the values derived from carbon, sodium, magnesium, and potassium — 
rather than much lower, as seen toward 23 Ori. In the main clouds toward C, Oph, however, most 
elements — even sulfur and calcium — give in the range 0.19-0.41 cm~^, and only iron gives 
a low value (0.03 cm~^). While including grain-assisted recombination (right panel) gives slightly 
better agreement for silicon, sulfur, manganese, and iron toward 23 Ori, the agreement among the 
various elements toward C, Oph is poorer (except that the value for iron is more consistent). In 
general, the values for the depleted elements (with s > 0) are increased, relative to those for the 
lightly depleted elements (with s ~ 0). Overall, the patterns of relative Ue seem rather similar for 
the two cases. 

If viewed by element instead of by sightline, the electron densities inferred from calcium are 
often higher than the values derived from other elements [but not toward C, Oph, unless grain- 
assisted recombination is important (for other elements)]. The values obtained from iron (and 
perhaps manganese and nickel) are often low (but not toward 1 Sco or 5 Sco); some apparent 
deficiencies of Mn I and Fe I had previously been reported by Joseph, Snow, &: Scab (1989) and 
Snow, Joseph, & Meyer (1986). The values derived from sulfur are also often low (but not toward 
C, Oph). Unfortunately, there are currently only a few sightlines with reliable high-resolution data 
for both neutral and singly ionized species for many elements. 

The differences in the element-to-element pattern of relative for these five lines of sight 
suggest that the differences in electron density inferred for any given line of sight are not just 
due to a single process or factor — but instead result from multiple additional processes that 

selectively affect the ionization balance of the various elements, to different degrees in the different 
lines of sight. The similarities between the patterns seen for (23 Ori, C, Ori), and also for (1 Sco, 
5 Sco) — and the differences between the Orion, Scorpius, and C, Oph patterns — are suggestive 
of regional variations in some (as yet undetermined) physical parameter (s). Fitzpatrick h Spitzer 
(1997), Welty et al. (1999b), and Weingartner &: Draine (2001) have all described some possible 
factors and/or processes - including problems with the adopted rates, dielectronic recombination, 
dissociative recombination, density-dependent depletions, charge exchange with small grains, and 
charge exchange with protons — but consistency among the various estimates for remains elusive. 

4-3.2. Dissociative Recombination of CH^ as a Source of C I 

In principle, dissociative recombination of CH+ "should" have a measurable impact on the 
ionization balance of carbon. The dissociative recombination rate has been estimated as ~ 3.3 
X 10-^ cm^ s-i at T ~ 100 K (Mitchell & McGowan 1978). Since A^(C Il)/iV(H) ~ 1.4 x 10"^ and 
iV(CH'*")/Ar(H) ranges from about 1-25 x 10~^ (from examination of various references), the ratio 
[aa Ar(CH+)]/[ar iV(C II)] — which gives the relative contributions of dissociative and radiative 
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recombination to C I — could be between 0.3 and 6.6. If dissociative recombination were generally 
dominant, we would expect that A^(C l)/iV(X I) would be proportional to A^(CH+)/iV(H), for some 
element X whose depletion is roughly constant. In Figure 8, we show the relationships between those 
ratios, for X = K, Na, and S. Of those ratios, consideration of the ionization potentials of the neutral 
species suggests that A^(C I)/A^(S I) should be the least susceptible to any effects of differences 
in the far-UV extinction. While the correlation probabilities (Table 4) indicate that there is no 
correlation for S I or K I, the two methods which include limits suggest a possible correlation for 
Na I. The regression fit, however, yields a slope consistent with no correlation (Table 5); including 
the available limits seems unlikely to increase it significantly. Probable differences in distribution 
between CH+ and hydrogen might act to obscure any correlations, however. As noted by Welty 
et al. (1999b) for 23 Ori, CH+ may not be entirely coextensive with the various trace neutral 
species, and so would not have as strong an impact on the ionization balance of carbon. Inclusion 
of grain-assisted recombination will reduce the relative importance of dissociative recombination 
and (probably) add scatter in comparing the ratios considered in this section — making it more 
difficult to discern the possible effects of dissociative recombination. 

Even if dissociative recombination of CH"'" is not the dominant source of C I, however, it may 
play a significant role in some clouds. Table 6 and Figure 6 both indicate that the electron density 
derived from carbon may be slightly higher than the values obtained from sodium and potassium 
for 23 Ori and ( Oph [where A^(CH+) is greater than 10^^ cm~^], but is slightly lower for ( Ori, 
1 Sco, and S Sco [where A/"(CH+) is < 3 x 10^^ cm"^]. [The shallow far-UV extinction toward S 
Sco and other stars in that region may also contribute to lower A^(C I), however (Welty 1989).] 
Enhanced C I toward some stars in the Pleiades (WH) may also be related to the relatively high 
iV(CH+) there (White 1984; Crane, Lambert, k Sheffer 1995). 

4-3.3. Does He Depend on nn? 

Different sources for the electrons in diffuse interstellar clouds have corresponding different de- 
pendences of Ue on the local hydrogen density nn. Those dependences, in turn, would yield different 
relationships between the local densities (and corresponding column densities) of various trace neu- 
tral species and those of hydrogen. Table 7, lists local hydrogen densities, electron densities, and 
fractional ionizations (ng/nn) estimated for the main clouds toward twelve stars; values for other 
discernible component groups are also given for several of the sightlincs. The estimates for nn are 
obtained from thermal pressures (uyiT) derived from analyses of the C I fine-structure excitation 
equilibrium (Jenkins & Shaya 1979). The C I column densities are derived from high-resolution 
(FWHM ~ 3 km s~^) UV spectra obtained with HST, in order to obtain accurate values for specific 
component groups (Wannier et al. 1999; Welty et al. 1999b and unpublished; Zsargo 2000; Jenk- 
ins & Tripp 2001; Sonnentrucker et al. 2002 and in prep.). The corresponding temperatures are 
obtained from analyses of the rotational excitation of H2 (Savage et al. 1977; Welty et al. 1999b; 
Rachford et al. 2002); values in parentheses are somewhat arbitrarily assumed (as 100 K). The last 



-25- 



four columns give the observed N{Ca l)/N{Ca II) ratio, then two sets of values for the electron den- 
sity inferred from that ratio, the average electron density derived from the corresponding ratios for 
carbon, sodium, and potassium (assuming typical depletions where necessary), and the fractional 
ionization using that average Ue- The first set assumes only radiative recombination; the second set 
includes both radiative and grain-assisted recombination, as described above (§§ 4.2.1 and 4.3.1). 
Electron densities in parentheses have been estimated from total sightline column densities — and 
so should be viewed as lower limits to the values in the gas traced by the neutral species. For all 
the electron density estimates, we use radiative recombination coefficients (a,.) corresponding to 
the adopted temperature. We assume the average (WJl) interstellar radiation field for all cases 
except the three Scorpius sightlines, where we multiply that field by a factor 5 (one possible expla- 
nation for the lower than usual abundances of some trace neutral species seen in those sightlines; 
see WH). The assumption of the unattenuated field should give at least rough estimates for Ue in 
these relatively diffuse clouds [with individual component E{B — V)< 0.2]. 

The values for nn, rig, and Ue/nji collected in Table 7 reveal several interesting trends. In 
the upper panel of Figure 7, we plot the values of He determined from calcium (squares) and 
from carbon, sodium, and potassium (CNaK) — with and without grain-assisted recombination 
(circles and triangles, respectively) — versus the local density nn- The open symbols denote the 
three sightlines in Scorpius, where the radiation field was assumed to be 5 x stronger. For this 
relatively small sample, the density ranges from 6-335 cm~^ (a factor greater than 50), but ne(rad) 
(CNaK) ranges only from 0.04-0.23 cm~^ (a factor of 6) — with no apparent dependence of ne(rad) 
on nn- The electron densities computed when grain-assisted recombination is included, however, 
appear to be somewhat smaller for larger nn — reflecting the increased effect of charge exchange 
at higher densities. On the other hand, the ne(rad) derived from calcium (which is not significantly 
affected by grain-assisted recombination, if the sticking factor s = 1) may show slight increases with 
increasing density — but this result depends somewhat on the assumed stronger radiation fields 
for the Scorpius sightlines. We note, however, that the ratio ne(Ca)/ne(CNaK) is also higher for 
those three sightlines (relatively independent of the strength of the radiation field). Similar trends 
are seen versus the thermal pressure nuT, as the temperature ranges only over a factor of 2 for this 
small sample. While the sample of component groups in Table 7 is likely not entirely representative 
of the diffuse ISM [a number of the groups have thermal pressures much higher than the typical 
values (Jenkins &; Tripp 2001)], these trends appear to be present even when the higher pressure 
groups are excluded. 

The apparent lack of dependence of ne(rad) on and the possible inverse relationship between 
He and nn when grain-assisted recombination is included are rather surprising. If the electrons are 
due primarily to photoionization of carbon and other abundant heavy elements, wc would expect 
Ue oc nn- If the electrons instead are due primarily to cosmic ray ionization of hydrogen, we would 

1 /2 

expect rig oc (e.g., Hobbs 1974; Kulkarni & Hcilcs 1987). From the discussions above (§§ 4.2.1 
and 4.3.1), we also would expect ne(rad) = rig + (ag nE_)/ar for s = (as adopted for carbon, 
sodium, and potassium), with the second term often dominant. We note, however, that some of the 
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models of Liszt (2003) — particularly those with enhanced cosmic ray ionization; mild depletion of 
carbon, nitrogen, and oxygen; and less than the nominal grain-assisted recombination — predict 
relatively constant Ug for nn ^ 100 cm~^. 

A lack of dependence of on rifj would also make it difficult to understand some of the 
observed relationships between the column densities of X I, X II, and H, for the typically severely 
depleted elements X = Fe and Ca. If radiative recombination dominates for iron and calcium 
(with s ~ 1), then we would expect that n{X'^) oc Ug n(X+) oc Ug nn ^(X). The relationships 
reported by Wakker & Mathis (2000) suggest that N{Fe II) <x [N{K)]^-^ and A^(Ca II) oc [A^(H)]0-2 
— presumably reflecting increasingly severe depletions in higher column density (denser?) gas. 
Combining column densities for Fe I and Ca I from this paper with A'^(H) from WH, we obtain 
iV(Fe I) oc [Af(H)]i-37±0-20 and N{Ca I) oc [iV(H)]i-25±o.i4^ Taking the ratios of the neutral and 
singly ionized species, we thus find that A^(X I)/A^(X II) — usually assumed to be proportional to 
Ug — is roughly proportional to A^(H) for both calcium and iron. 

Hobbs (1974) had noted that the observed roughly quadratic dependences of A''(K I) and 
A^(Na I) on A^(H) would be consistent with the essentially constant fractional ionization expected 
if most electrons came from photoionization of carbon, and concluded that cosmic ray ionization of 
hydrogen [which would yield shallower relationships between those N{X I) and A''(H) under the as- 
sumption of photoionization equilibrium] must not be the main source of electrons. If grain-assisted 
recombination dominates the production of the trace neutral species for very mildly depleted ele- 
ments like sodium and potassium (with s ~ 0), however, then we would expect that n(X.^) oc nn 
n(X+) oc 5(X) — consistent with the observed roughly quadratic relationships. Those observed 
quadratic relationships thus would not exclude significant cosmic ray ionization of hydrogen in 
diffuse interstellar clouds (Liszt 2003). 

4.3.4- Fractional Ionization 

The other puzzling aspect of the ionization balance in the main clouds toward 23 Ori is the 
generally high fractional ionization — of order 0.01 — implied by the adopted Ug and the nn 
obtained from analysis of the C I fine structure excitation. That fractional ionization is a factor 
50 higher than the value ~ 2 x 10~^ expected from photoionization of carbon and other abundant 
heavy elements, and is still a factor of order 10 higher than the values obtained by including 
ionization of hydrogen and helium due to cosmic rays and x-rays, even if that ionization rate is 
as high as 10~^^ s~^ (e.g., Kulkarni & Hciles 1987; Wolfire et al. 1995; Wcingartncr & Draine 
2001). Welty et al. (1999b) had conjectured that charge exchange with large molecules might have 
enhanced the trace neutral species by perhaps an order of magnitude (based on the models of Lepp 
et al. 1988), but the more detailed calculations of Weingartner &: Draine (2001) suggest that the 
effects of charge exchange will be somewhat smaller (at least in the main clouds toward 23 Ori; 
see Table 6). In addition, Weingartner &; Draine note that grain-assisted recombination will also 
affect the abundance of H"^ in predominantly neutral regions — countering the effects of cosmic 
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ray and x-ray ionization of hydrogen to some extent (see also Liszt 2003). The question (again) is 
how representative the clouds toward 23 Ori are of the diffuse ISM. 

In the lower panel of Figure 7, we plot the fractional ionizations for the component groups listed 
in Table 7 versus nn, using the same symbols as in the upper panel. While ng/nn is often assumed 
to be constant, the values for the small sample shown in the figure appear to decrease at higher 
densities (corresponding to the trends seen for rig in the upper panel). Nine of the 13 component 
groups have ne(rad)/nH > 10~^ (triangles); three of the groups (those toward HD 185418, ( Per, 
and 23 Ori) exhibit values greater than 10~^. Inclusion of grain-assisted recombination (circles) 
reduces the values for three of those nine component groups to within the range 2-10 x 10~^ 
"expected" if both photoionization of carbon (and other abundant heavy elements) and cosmic ray 
ionization of hydrogen are significant sources of electrons, but also reduces the other four values 
originally within that range to well below 2 x 10~^. 

4.3.5. Electron Densities from N (Ca II) /N (Ti II) 

In connection with a large survey of interstellar Ti II absorption. Stokes (1978) developed 
a method for estimating rig from the ratio N{Ca ll)/N{Ti II) (see also Stokes & Hobbs 1976). 
This method would be of particular utility for the thinner diffuse clouds, where Ca II is likely a 
trace species (with Ca III dominant) and where many of the various trace neutral species are more 
difficult to detect. The method assumes both photoionization equilibrium and a constant ratio r = 
[^0(Ca)(^(Ca)]/[A0 (Ti)5(Ti)] in the interstellar gas [as suggested by the good correlation observed 
between A^(Ca II) and iV(Ti II)]. The value of r is calibrated by the N{Ca l)/iV(Ca II) ratio (or the 
corresponding rig) where Ca I has been detected. We then have ng ~ (T/ar) {r[N{Ti Il)/A^(Ca II)] 
— 1}~^, where the T/ar is for the Ca Il-Ca III equilibrium. 

The above discussions suggest several potential difficulties with this method, however. First, 
the assumption of photoionization equilibrium may not hold, as other processes may be significant 
(though grain-assisted recombination may not be important for these typically severely depleted 
elements). Second, calibrating r with the from calcium will yield a value that is too small (and 
thus estimated electron densities that are too large). Finally, in order to estimate rig (given r), one 
must know the temperature (in order to obtain the recombination rate Ur for Ca"*"^ + ng). If the 
temperature can be as high as several thousand K for warm, neutral gas (where much of the Ti II 
may reside), then could be as much as an order of magnitude smaller than the value for T ~ 
100 K which is usually assumed. 

4.4. Density Dependence of Depletions 

As noted above, trends in average linc-of-sight depletions with average density < riH > and with 
/(H2) have suggested that the depletions may depend (to varying degrees) on the local hydrogen 
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density nn- Using column density data collected from the literature and/or derived from new 
observations of stars in the Cepheus 0B3 association, Cardelli et al. (1991) found an inverse 
relationship between the ratios A^(CN)/A^(CH) and A^(Ca l)/A^(CH). Cardelli et al. surmised that 
this somewhat unexpected anti-correlation between Ca I and CN might be due to increasingly 
severe depletion of calcium in the denser gas responsible for the CN. Using a simple model for 
the gas-phase chemistry and assuming that photoionization equilibrium holds for calcium and that 
most electrons come from ionization of carbon, they concluded that the calcium depletion varies 
roughly as n^"^. This method thus seems able to provide specific relationships between the local 
density and the depletions of various elements X, if column densities for X I, CH, and CN are 
available for a sufficient sample of sightlines. Again, however, the above discussions suggest that 
the assumptions of photoionization equilibrium and of carbon being the primary source of electrons 
may not be valid. 

Using the larger sample of column densities now available, we may examine such relationships 
for several different trace neutral species (Figure 9). The correlation probabilities obtained from the 
survival analysis routines suggest that (anti)correlations may be present for Ca I and Fe I (Table 4) 
— though the values obtained from the detected points alone suggest otherwise. The relationship 
between iV(CN)/iV(CH) and iV(Ca l)/iV(CH) seems both steeper (slope < -2.5 rather than -1; 
Table 5) and less tight than the relationship plotted by Cardelli et al. (1991). Within the Cardelli 
et al. framework, steeper relationships correspond to weaker dependences of the depletion on nn; a 
slope of —2.5 would imply that the calcium depletion is proportional to n^^'^. Substituting Fe I for 
Ca I also yields a steeper slope (< —2.0), with a slightly weaker (anti) correlation than for Ca I. For 
A^(CN)/A^(CH) vs. A^(K I)/A^(CH), we find an essentially infinite slope [omitting the points for 
HD 37903 and HD 73882 with anomalously small A^(K I)/A^(CH)], since the latter ratio appears to 
be roughly constant for most lines of sight (WH). The analog of equation 10 in Cardelli et al. for 
potassium would thus imply that the depletion of potassium is proportional to n^^. The roughly 
constant ratio N{C T)/N{K I) (WH) and the apparently uniform depletion of carbon over wide 
ranges in < nn > and /(H2), however, suggest that the depletion of potassium does not depend on 
density, at least for the range in density probed so far.^^ 

Another way of gauging the possible dependence of depletions on density is to compare ratios 
such as A^(Ca I)/A^(K I) [oc (5(Ca)] with the ratio A^(CN)/A^(CH). Because relatively dense gas is 
necessary for the formation of CN, the ratio A^(CN)/A^(CH) should be a density indicator; simple 
chemical models suggest that n(CN)/n(CH) is proportional to (Cardelli et al. 1991). The plots 



^"While Cardelli et al. (1991) considered upper limits essentially as detections, we have chosen to fit only the 

points corresponding to detections - - which may contribute to the steeper slope derived in this paper. For all three 
trace neutral species considered here, a slope of infinity yields a better fit to the available data than the slope of — 1 
found for Ca I by Cardelli et al. 

^^The possible upturn in A'^(C I) vs. N{K I) for the highest column densities in the current sample may reflect 
uncertainties in the /-values for some of the C I lines (Jenkins & Tripp 2001) and/or the effects of UV extinction 
(Jenkins & Shaya 1979; Welty 1989). 
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of iV(Ca I)/N{K I) and N{Fe l)/N{K I) vs. N{CN)/N{CR) shown in Figure 10 reveal (at best) 
relatively weak anti-correlations (Table 4); only the correlation tests which include limits suggest 
that correlations may be present. The slopes computed for the detected points are marginally less 
than zero (Table 5); the available limits are consistent with those fits. Similar behavior is observed 
if Na I is substituted for K I (but for fewer lines of sight). If n(CN)/n(CH) goes as and if the 
depletions of sodium and potassium do not depend on the density, then these weak anti-correlations 
would be consistent with more severe depletion of calcium and iron in denser gas — but with a 
much weaker dependence on density (at most as n^^'^) than that found for calcium by Cardelli et 
al. 

Finally, for the small sample of component groups listed in Table 7, we may directly compare 
iV(Ca Il)/iV(H) with rin. In ah cases, the iV(Ca l)/N{Ca II) ratios suggest that Ca II is the 
dominant form of calcium in those component groups, so that A^(Ca II)/A^(H) oc 6(Ca). A weighted 
fit to the 11 data points yields the relation 6{Ca,) oc [nH]"^"^"'^; any differences in distribution 
between Ca II and H would likely make the dependence somewhat shallower. 

4.5. Individual Lines of Sight 

4.5.1. Cepheus OB 2 

Many of the strongest Ca I lines observed in this survey are toward stars in the Cepheus 0B2 
association. In several of the Cep 0B2 sightlines — toward HD 207198, u Cep, and A Cep — there 
is a strong, narrow Ca I component, with relatively high A^(Ca I)/A^(K I), at or near the red- ward 
edge of the absorption (though at somewhat different velocities in the three cases). The red-most 
Ca I component toward A Cep, at -y —2.7 km s"-^, also has one of the highest A'"(Ca I)/A^(Ca II) 
ratios in the sample. Analysis of the C I fine-structure excitation toward A Cep suggests that the 

— 2.7 km s""*^ component is characterized by a somewhat greater thermal pressure than the other 
strong, low-velocity components, which have log nuT (cm~'^K) of order 3.3 (Table 7; Jenkins &: 
Tripp 2001). 

Jenkins & Tripp (2001) find even higher pressures for the component near —35 km toward 
A Cep. The C I analysis yields pressures ranging from log n-nT ~ 5.3 for T = 160 K to log uhT 
~ 5.8 for T = 80 K; similar analysis of the O I fine-structure excitation yields log ~ 4.8. 
The relatively narrow width of that component as seen in K I — 6 ~ 0.8 km s"^ (WH) — implies 
that the temperature must be less than about 1500 K; if the turbulence in the gas is just sonic, 
then T would be of order 150 K. Intriguingly, the A^(Ca l)/A^(Ca II) ratio — and thus rig/ (T /or)'^ 

— in the —35 km s~^ component is not much greater than the values found for the low- velocity 
components — even though the pressure is at least 30 times higher. 
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4-5.2. Orion 

The narrow blue-most Ca I components toward l Ori, e Ori, and ^ Ori (at slightly different 
velocities) are quite strong, relative to the corresponding narrow components in Na I and K I 
(WHK: WH; Price et al. 2001). With Na I 6- values of order 0.35-0.45 km s~^, these components 
must have T < 275 K. Toward e Ori, the component at 3.1 km s~^ has a high N{Ca l)/N{Cat II) 
ratio — nearly 8 times higher than the value for the main components near 24 km s~^. Toward 
Ori, however, the component at —1.8 km has a somewhat lower A^(Ca l)/A^(Ca II) ratio, 
more comparable to that in the main components near 24 km s^^. We note, however, that this 
component is apparently time- variable (Price et al. 2001), and that the corresponding A^(Ca II) may 
actually be significantly smaller than the value used here, as the Ca II component at that velocity 
(in the published profile fits) has a much higher 6- value (WMH; Price et al. 2001). Analysis of the 
C I fine-structure excitation toward (" Ori indicates that the thermal pressure in the —1.8 km s~^ 
component is more than an order of magnitude larger than that in the main components (Tabic 7). 
Jenkins & Peimbert (1996) found a systematic broadening, with increasing rotational level J, of 
the H2 lines near v ~ km in IMAPS spectra of ( Ori, and proposed that the H2 is formed 
in cooling, compressed post-shock gas (see also Welty et al. 2002). The gas-phase abundances 
(relative to S and Zn) in the gas near v = km s"^ are very similar to those typically found in 
warm, quiescent, diffuse clouds, however — i.e., there are no obvious enhancements of typically 
severely depleted elements that would suggest that that gas has been shocked. 

4.5.3. SN 1987A 

Magain (1987) found unusually strong Ca I absorption in a number of velocity components 
toward the LMC SN 1987A. In particular, the ratios 7V(Ca l)/iV(Na I) and N{C& l)/N{Mg I) 
in some of the component groups are significantly higher than the values commonly found in the 
Galactic ISM; the A''(Ca I)/A/"(Ca II) ratio, however, generally is closer to the Galactic values 
(Welty et al. 1999a). While much of the apparent enhancement of Ca I toward SN 1987A is likely 
a result of less severe calcium depletion, dielectronic recombination of Ca II in relatively warm gas 
(r ~ 4500 K) may also be a factor for some of the intermediate-velocity components (Welty et al. 
1999a). 

5. Summary 

We have analyzed high-resolution (FWHM ~ 0.3-1.5 km s^^) spectra of interstellar Ca I A4226 
absorption toward 30 Galactic stars. Typical 2a equivalent width limits for narrow individual 
components are 0.16-0.95 mA, which correspond to Ca I column density limits of 0.6-3.6 x 10^ 
cm~^. These spectra are thus characterized by higher resolution and sensitivity than the Ca I 
spectra previously available for most of these lines of sight. 
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Voigt-profile fits to the observed fine profiles yielded a total of 112 individual Ca I components 
— which generally correspond to components seen in Ca II, Na I, and K I. The median Ca I 6- value 
for the highest resolution spectra, 0.66 km s^^, is comparable to those found for Na I (WHK) and 
K I (WH), as would be expected if the three species are essentially coextensive. The generally 
larger 6-values found for Ca II (WMH) suggest that Ca II may be more broadly distributed, in 
gas characterized by somewhat higher temperatures and/or turbulent velocities. For corresponding 
individual components at the same velocity, the ratio A^(Ca I)/A^(K I), indicative of the calcium 
depletion, ranges over a factor of about 500; the ratio A^(Ca l)/A^(Ca II), equal to ne/(r/a,.) in 
photoionization equilibrium, ranges over a factor of about 50. Comparisons between Ca I and 
various other species yield the following results: 

• If photoionization and radiative recombination dominate the calcium ionization balance, then 
the observed column densities of Ca I and Ca II suggest that Ca II is the dominant form of 
calcium wherever Ca I is detected. Ca III could still be dominant, where A'^(Ca l)/N{Ca II) is 
smallest, if charge exchange with small grains is more important than radiative recombination. 

• The ratio N{Ca l)/N{Ca II) does not seem to depend on the fraction of hydrogen in molecular 
form /(H2) — suggesting either that ne/{T/ar) does not depend on tzr or that /(H2) is not 
a reliable indicator of density or that photoionization equilibrium does not apply. 

• Ca I seems to be systematically more abundant, relative to other trace neutral species, than 
would be expected if photoionization and radiative recombination control the ionization bal- 
ance of heavy elements. The typically small Ca I 6- values imply that the temperature cannot 
be high enough (> 3000 K), in general, for dielectronic recombination of Ca+ to account for 
that apparent overabundance. Grain-assisted recombination is unlikely to significantly affect 
the ionization balance of the typically severely depleted calcium. 

• The electron densities inferred from the ratios A^(X l)/A^(X II), for different elements X, 
can differ by more than an order of magnitude along a given line of sight, if photoionization 
equilibrium is assumed. In particular, the Ue obtained from N{Ca l)/N{Ca II) is often 
significantly higher than the values derived from other elements. The values obtained from 
the corresponding ratios for sulfur, manganese, iron, and nickel are often lower than those 
derived from carbon, sodium, and potassium. Those differences are not eliminated if charge 
exchange with small grains dominates the production of (at least) some of the trace neutral 
species X I. The patterns of relative Ue for different lines of sight show both similarities 
and differences — suggesting that multiple additional processes besides photoionization and 
radiative recombination significantly affect the ionization balance of heavy elements — to 
different degrees in different lines of sight. 

• The lack of a strong correlation between N{C l)/N{X I) and A^(CH+)/Ar(H), for X = Na, K, 
and S, suggests that dissociative recombination of CH^ is generally not the dominant source 
of C I — though it may make a significant contribution in some individual cases. 
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• For a small number of component groups with independent estimates for Ue and ran, the 
resulting fractional ionization is often much greater than the value 2 x 10~^ that would 
be expected if most of the electrons come from photoionization of carbon and other heavy 
elements — suggesting that rie has been overestimated (if only radiative recombination is 
assumed) and/or that hydrogen may be slightly ionized in predominantly neutral clouds. 
Inclusion of grain-assisted recombination can reduce the apparent discrepancies in some lines 
of sight. 

• For that same small sample, the average electron density for carbon, sodium, and potassium 
(assuming only radiative recombination) appears to be roughly constant — at ng ~ 0.14±0.07 

cm~^ (mean it standard deviation) independent of the local hydrogen density nu- If grain- 
assisted recombination is included, the resulting rie appear to be smaller for larger njj — 
reflecting the increased effects of charge exchange at higher densities. Observed relationships 
between the column densities of various neutral and singly ionized species and A^(H) are 
difficult to understand if ng is either independent of or inversely proportional to nn- 

• Relationships between ratios of the column densities of Ca I, K I, Na I, CN, and CH suggest 
that the depletion of calcium has a much weaker dependence on density than that inferred 
from previous comparisons based on more limited data. 



It has been common practice to assume photoionization equilibrium to estimate the electron 

density in diffuse interstellar clouds, then to assume a fractional ionization equal to the gas-phase 
carbon abundance to estimate nn- The results discussed above indicate that such estimates must 
be viewed with caution — especially if they are based on only one ratio iV(X l)/A^(X II). Because 
lines from both Ca I and Ca II are available in the optical, the Ca ionization balance is often used — 
but the resulting Ue can be higher than the values inferred from other elements by nearly an order 
of magnitude. (And the "true" values from calcium would be even higher if Ca II is more widely 
distributed than Ca I.) Even if other elements are used to determine ng, the resulting fractional 
ionization (using nn determined from C I fine-structure excitation) can still be significantly larger 
than the gas-phase abundance of carbon if only radiative recombination is considered. Inclusion 
of grain-assisted recombination appears to yield more consistent results for Ue and ne/riu in some 
cases, but not for every element or for every line of sight. It may prove interesting to revisit some 
of the possible correlations explored in this paper when more detailed and extensive abundances for 
the various atomic and molecular species become available for more individual interstellar clouds. 
High-resolution observations of lines from many different ncutral/first-ion pairs, for additional 
sightlines sampling different environmental conditions, may help us to understand the various 
processes affecting the ionization balance in diffuse clouds. 
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A. High- Resolution Spectra of Fe I 

During our 1997 March AAT run, we also obtained high-resolution spectra (FWHM ~ 0.32 
km s~^) with the UHRF of the weak Fe I line at 3719A for nine lines of sight; Barlow et al. (1995) 
obtained similar spectra for (" Oph. The instrumental setup, data reductions, and analysis were 
very similar to those described above for the Ca I spectra. The normalized Fe I spectra for six of 
the sightlines are shown in Figure 11, together with the profiles of either Ca I or K 1. The measured 
Fe I equivalent widths and the component structures derived from fits to the line profiles are given 
in Table 8; equivalent widths from the few previous studies to have measured this line are given in 
a footnote to the table. Wc have used / = 0.0413 (Morton 1991) in deriving the column densities. 
Examination of the profiles in Figure 11 indicates that the component structures seen in Fe I arc 
very similar to those for Ca I; the Fc I 6-valucs arc similar to those found for the corresponding 
components in Ca I and K I. There are obvious variations in the ratio A^(Fe l)/A^(Ca I), however 
(compare ^ Ori with S Sco, a Sco, and p Oph) — reminiscent of the variations in relative derived 
from iron and calcium in different sightlines (§ 4.3.1 and Fig. 6). 

In Figure 12, we plot A'^(Fe 1) vs. N{Ca. 1) for all lines of sight with data for both species 
(see below). The slope of the best-fit (solid) line is about 1.13±0.14 — marginally steeper than 

that for a linear relationship between the two species, and perhaps suggestive of slightly less severe 
depletion of iron, relative to that of calcium, in the higher column density sightlines. The scatter 
is about 0.21 dex, similar to that for A'^(Ca 1) vs. A^(K I). 

B. Column Densities for Ca I, Ca II, Fe I, and S I 

In Table 9, we list total line-of-sight column densities for Ca I, Ca II, Fe I, and/or S I for 135 
sightlines for which data for Ca I are available. An asterisk in the third column denotes sightlines 
included in the appendix to WH, where similar data for H, H2, Na I, K I, Li I, C I, and CH are 
collected. We have used /-values from Morton (1991), with updates for some lines as listed by 
Welty et al. (1999b), for determining the column densities. 

Ca I: The Ca I column densities for sightlines observed for this paper are based on the fits to 
the high- resolution line profiles. Values for some more heavily reddened sightlines are from a high- 
resolution survey of optical absorption lines toward stars in the FUSE PI team program studying 
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H2 in translucent cloud sightlines (Welty et al., in prep.; see Rachford et al. 2002). Values for other 
sightlines were derived from equivalent widths of the A4226 line (see references); thirty of those 
were measured from high S/N spectra obtained with the ARC cchcllc spectrograph (FWHM ~ 8 
km s^^) for a survey of diffuse interstellar bands (J. Thorburn, priv. comm. ). For Wa < 15 mA, 
we have assumed that the lines are optically thin, and have given lower limits for stronger lines. 

Ca II: We have only included Ca II column densities derived from high-resolution spectra of 
the A3933 K line (FWHM < 2 km s"^; mostly from WMH), as the Ca II lines can be strong and 
complex where Ca I is detected. 

Fe I: For Fe I, very few high-resolution optical spectra are available (this paper; Barlow et al. 
1995). We include values obtained from lower resolution spectra (e.g., HST/GURS echelle spectra 
of the lines at 2484 and 2523 A), as the measured absorption lines are generally relatively weak 
(e.g., Welty et al. 1999b). Forty-one of the values were derived from equivalent widths of the Fe I 
lines at 3719 and 3859 A, obtained from the ARC diffuse band survey spectra noted above (J. 
Thorburn, priv. comm.). 

S I: For S I, few results from high-resolution UV spectra have been reported, and there are no 
lines at optical wavelengths. We have used component structures derived from higher resolution 
K I and/or Na I spectra (when available) to model the measured S I equivalent widths (typically 
for the lines at 1295, 1425, and 1807 A; e.g., Welty et al. 1999b); otherwise, we have scaled column 
densities in the original references to the current /-values. 

While we have not listed uncertainties for the column densities in Table 9 (or in Table 7 
of WH), where possible we have tabulated, derived, or estimated uncertainties for those column 
densities for use in the correlation and regression analyses^^. For column densities derived from 
fairly weak absorption lines (Wa ^15 mA), the column density uncertainties generally correspond 
to the uncertainties in the equivalent widths — which are often of order 10-20% (0.04-0.08 dex) 
even for confidently detected lines. Most values for Ca I, Fe I, and Li I — and some values for 
Na I, K I, Ca II, and CH — are based on such weak lines. For column densities derived from 
very strong lines (many of the values for Na I, K I, and Ca II), the formal uncertainties in the 
total column densities derived from fits to high-resolution, high-S/N spectra are often < 0.02 dex. 
The true uncertainties are likely to be somewhat larger, however, as the fits assume the number 
of components present and (in many cases) fix some of the 6-values to facilitate convergence. 
Moreover, component structure information from other, similarly distributed species also may be 
used to fit very strong lines and/or absorption-line profiles observed at slightly lower resolution. 
While the effects of those assumptions are difficult to quantify, in some cases the column densities 
can be compared with those derived from weak lines of the same species. For example, K I column 
densities derived from fits to the strong A7698 line were within 0.05 dex of those obtained from the 
weak A4044 line for three of the five sightlines with reliable data for both lines (WH) . We therefore 
have adopted a minimum uncertainty of 0.03-0.04 dex for column densities derived from fits to 
relatively strong lines (for all species). For the species considered in this paper, the typical and 
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adopted minimum uncertainties are: C I (0.2/0.1 dex), Na I (0.04/0.04 dex), S I (0.1/0.05 dex), 
K I (0.03/0.03 dex), Ca I (0.06/0.03 dex), Ca II (0.03/0.03 dex), Fe I (0.07/0.04 dex), H I (0.1/0.05 
dex), H2 (0.2/0.05 dex), CH (0.04/0.03 dex), CH+ (0.05/0.03 dex), and CN (0.09/0.04 dex). 



^^In some cases, the original references (including our previous surveys) do not give uncertainties for either the 
measured equivalent widths or the derived column densities. In other cases, it is not clear what effects were included 
(e.g., photon noise, continuum placement, background uncertainties), or exactly how the listed uncertainties were 
calculated — making it difficult to compare values quoted in different studies. An updated compilation of column 
densities (with estimated uncertainties) may be found at http://astro.uchicago.edu/home/web/welty/coldens.html. 
Compilations of equivalent widths measured from high-resolution optical spectra (with uncertainties) may be 
found at http://astro.uchicago.edu/home/web/welty/ew-atom.html (for Na I, K I, Ca I, and Ca II) and 
http://astro.uchicago.edu/home/web/welty/ew-mol.html (for CN A3875, CH A4300, and CH+ A4232). 
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Fig. la. — The interstellar Ca I, K I (or Na I), and Ca II profiles. The vertical scale has been 
expanded in some cases to show weaker lines (especially of Ca I) more clearly (noted by an asterisk 
near both vertical axes). The sources of the spectra have been noted at the right, just above the 
continuum [M = McDonald 2.7m coude (FWHM = 0.53-0.56 km s'^); K = KPNO coude feed 
(FWHM = 1.2-1.5 km s"!); A = AAT UHRF (FWHM = 0.3-0.4 km s'^); E = ESO 3.6m + CES 
(FWHM = 1.2 km s'^); H = Hobbs (PEPSIOS; FWHM = 1-2 km s'^)]. Tick marks above the 
spectra indicate the components found in fitting the profiles. The solid triangles denote vlsr = 
km s""-*^. 
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Fig. 2. — Fits to the Ca I profile observed toward o Per. The top panel shows a two-component fit 
— which is not adequate; the bottom panel shows the adopted four-component fit. In both panels, 
the individual components (with velocities given by the tick marks) are shown by dotted lines; the 
composite fits are shown by the smooth solid line; the data are shown by the upper solid histogram; 
the residuals (data minus fit) are shown by the lower solid histogram. 
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Fig. 3. — N(Ca I) vs. N(K I). Open circles denote Sco-Oph sightlines; open triangles denote other 
"discrepant" sightlines identified by WH. The size of each plus sign indicates the iblcr uncertainties 
in each quantity. The slope of the best-fit line (solid) is 0.60ib0.06 — significantly shallower than 
that for a linear relationship (dotted line) — which may reflect increasing severity of Ca depletion 
in higher column density (higher mean density?) sightlines. 
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Fig. 4.— N(Ca I)/N(K I) and N(Ca I)/N(Ca II) vs. /(Hg). Open circles denote Sco-Oph sightlines; 
open squares denote Orion sightlines; filled circles denote other sightlines. Crosses (connected by 
lines) denote values for (some) individual components, assumed to have the same /(H2) as the 
average for the line of sight. The total line-of-sight ratios probably underestimate the true values 
where Ca I is detected over a much smaller velocity range than K I and/or Ca II. N(Ca I)/N(K I) is 
smaller for larger /(H2) — consistent with the ratio being an indicator for D(Ca). N(Ca I) /N(Ca II) 
shows no obvious dependence on /(H2) — so either ne/{T/a) is independent of nn or /(H2) is not 
a reliable indicator of nn or photoionization equilibrium does not hold. Typical la uncertainties 
in both ratios are of order 0.07 dex. 
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Fig. 5. — N(Ca I)/N(X I) vs. N(Ca II)/N(H) for various elements X. Open circles denote Sco- 
Opli sightlines; open squares denote Orion sightlines; open triangles denote main LMC components 
toward SN 1987A. The size of each plus sign indicates the ±1(7 uncertainties in each quantity. The 
solid lines give the expected (linear) relationships if X is undepleted (assuming the photoionization 
and radiative recombination rates in Tabic 3); the dotted lines give the relationships for typical 
depletions (—0.4 dex for C, —0.6 for Na, —0.7 for K, 0.0 for S); for comparison, the dashed lines 
show the fits to the data with the slopes fixed at unity. 
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Fig. 6.— {left) Logarithm of Ue = N(X I)/N(X II) * (r/a^), for T = 100 K and the WJl radiation 
field, relative to no (the average value for C, Na, and K), for the main components toward five 
stars (Table 6). Open circles denote cases where X II was estimated based on N(H) and typical 
depletions; representative 1 a error bars are shown for 23 Ori. Note the similarities between 23 Ori 
and ( Ori, and between 1 Sco and 6 Sco — and the differences between those two pairs and Oph. 
The points for S, Mn, Fe, and Ni are often low; Ca is often high, (right) The same plot, except that 
charge exchange with small grains is included (Weingartner Sz Draine 2001). The circles denote 
the Ue/riQ values for the adopted set of sticking parameters s (see text); the attached "error bar" 
denotes the range between all s=0 and all s=l. The assumed depletions of Na and K are less severe 
by ~ 0.2 dex. The element-to-element differences are generally similar to the case where radiative 
recombination is dominant, however. 
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Fig. 7. — Electron density (top) and fractional ionization ne/riu (bottom) versus local hydrogen 
density nn- Squares are values derived from calcium; circles and triangles are average values 
derived from carbon, sodium, and potassium — with and without grain-assisted recombination, 
respectively (see Table 7). Open symbols denote the three Scorpius sightlines, where the radiation 
field was assumed stronger by a factor of 5. Typical uncertainties are less than about 0.3 dcx for 
log(nH), but may be somewhat larger for log(ne). Note that ne(rad) (CNaK) (triangles) appears 
to be independent of riu, within those uncertainties. 
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Fig. 8.— N{C l)/N{X I) vs. A^(CH+)/iV(H), for X = Na, K, and S. The size of each plus 
sign indicates the iblfj uncertainties in each quantity. Dotted hues show the mean values of the 
A^(C I)/A^(X I) ratios; solid lines show the correlation fits, with slopes 0.0-0.2 (Table 5). Linear 
correlations (with slope = 1.0) would have been expected if dissociative recombination of CH+ 
dominates the production of C I. 
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Fig. 9.— iV(CN)/iV(CH) vs. iV(X I)/Af(CH), for X = Ca, K, and Fe (following Cardelli et al. 
1991). The size of each plus sign indicates the ±1(7 uncertainties in each quantity. For Ca and Fe, 
there may be inverse relationships between the two ratios, but with steeper slopes (< —2.5 and < 
—2.0; dotted lines) than the —1 found by Cardelli et al. for Ca — suggesting that the depletions 
of Ca and Fe do not depend as strongly on the local density as n^'^. For K, the slope is essentially 
infinite, as A^(K I)/iV(CH) is roughly constant for most Galactic sightlines (WH). 
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Fig. 10.— iV(Ca I)/iV(K I) and iV(Fe I)/N{K I) vs. iV(CN)/iV(CH). The size of each plus sign 
indicates the itlcr uncertainties in each quantity. Weak trends, with slightly negative slopes, may be 
present (especially if the limits are considered). If n(CN)/n(CH) ~ iV(CN)/A?^(CH) is proportional 
to n^, and if the depletion of K docs not depend on density, then those weak trends would suggest 
that the depletions of Ca and Fe go (at most) roughly as n^^'^. 
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Fig. 11.— The interstellar Fe I A3719 profiles obtained with the AAT/UHRF (FWHM ~ 0.3 
km s~^) toward six stars. The corresponding Ca I or K I profiles are given for comparison. Tick 
marks indicate components found in fitting the profiles. The vertical scales are the same for Fe I 
and Ca I in all cases; note the variations in relative strength of the two lines. 
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Fig. 12. — N(Fe I) vs. N(Ca I). The size of each plus sign indicates the itlo" uncertainties in each 
quantity. Open circles denote Sco-Oph sightlines. The slope of the best- fit line (solid) is 1.13±0.14 
— consistent with a linear (or perhaps slightly steeper) relationship between the two species (dotted 
line) . 
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Table 1. 


Ca I Equivalent Widths^ 


Star 





CF - - - 




IViCU - - - 


- - - u nrvr - - - 


• b 

Previous 




S/N 


Wa 
(mA) 


0/ IN 


Wa 
(mA; 


C /1\T W , 
0/ IN * • A 

(mA; 




K Cas 


205 


3.2±1.3 


155 


2.5±1.7 




<0.9 (1) 


HD 21291 


135 


5.7±1.4 








10.3±1.5 (2) 


o Per 






160 


2.0±0.4 




1.5±0.9 (3); <2.0 (4) 


C Per 






190 


1.6±0.4 




<0.9 (1); 0.9±0.6 (3); <2.4 (4) 


(, Per 


220 


1.6±0.5 


165 


1.5±0.8 




<1.1 (1); 2.3±0.6 (2) 


a Cam 


120 


1.9±1.0 








2.4±0.4 (5) 


23 Ori 


65 


4.1±1.3 








4.7±0.9 (2); 5.6±0.8 (5); 3.2±1.6 (6) 


\ Ori 


120 


2.1±0.8 


190 


2.9±0.6 




3.8±0.6 (5) 


t Ori 


185 


0.4±0.3 








0.7±0.2 (5) 


e Ori 










335 1.26±0.30 




C Ori 


145 


1.5±0.8 


250 


0.5±0.3'= 




1.0±0.2 (5) 


K Ori 


300 


0.3±0.2 










Ori 


150 


11.0±1.2 










HD 72127 A 










75 0.66±0.35 




a Vir 


335 


<0.2 










5 Sco 










225 0.66±0.29 


0.7±0.4 (3); 0.9±0.2 (5); 0.65±0.2 (6); <1 (7) 


Sco 










215 1.26±0.24 


<0.9 (1); 1.2±0.6 (3); 1.2±0.2 (6) 


fj Sco 










210 0.95±0.29 


1.0±0.4 (3); 1.1±0.4 (6); l.OiO.8 (7) 


p Opli 


80 


5.7±1.2 








4.8±0.8 (2); 4.5±1.2 (6); 5±2 (7) 


X Opli 


80 


2.1±1.3 








<0.6 (6) 


COph 


265 


0.8±0.4 


295 


0.8±0.2 




<1.5 (1); 1.3±0.8 (3); 1.6±0.8 (6); 1.3 (8) 


HD154368 


65 


5.3±3.8 








6.5±4.6 (9) 


P Cyg 


85 


9.3±2.8 










a Cyg 


255 


l.OiO.6 


325 


1.0±0.5 




<1.3 (1) 


HD207198 


55 


12.7±3.4 








15.3±1.3 (2); 15 (10) 


9 Cep 


95 


7.5±1.9 








12.1±1.9 (2); 12 (10) 


19 Cep 


65 


7.0±3.1 








9.3±2.7 (2); 7.2±1.0 (4) 


A Cep 






115 


10.4±1.8 




9.2 (1); 9 (in) 


V Cep 


80 


14.0±2.6 


125 


16.2±1.4 




<20.9±2.3 (2) 


6 Cas 


35 


4.3±2.6 











''Uncertainties and limits are 2a. 

''Equivalent width (reference) for previous data. 

'^Not including component at —2 km s~^. 



References. — 1 = Hobbs 1971; 2 = Welty (unpubl. CF camera 5); 3 = White 1973; 4 = Allen 1994; 5 = Federman & 
Hobbs 1983; 6 = Lambert & Danks 1986; 7 = Snow et al. 1986; 8 = Shulman et al. 1974; 9 = Gredel et al. 1993; 10 = 
Chaffee & Dunham 1979 



Table 2. Ca I Component Structures 



Star Obs^ Comp vq N 

(km s-l) (10'' cm-2) (km s"!) 



K Cas 


M 


1 


-25.07 


1.3 


1.27 






2 


-21.06 


4.6 


0.80 






3 


-13.86 


1.4 


0.91 






4 


-11.19 


0.8 


(1.00) 






5 


—6.44 


0.9 


1.02 


HD 21291 


K 


1 


-10.36 


6.1 


1.86 






2 


-6.20 


9.3 


1.19 






3 


-0.41 


6.1 


1.67 


o Per 


M 


1 


10.53 


1.0 


[0.35] 






2 


11.55 


1.8 


0.33 






3 


13.45 


3.8 


0.75 






4 


14.85 


0.7 


0.37 


C Per 


M 


1 


11.35 


0.7 


0.60 






2 


13.35 


2.9 


0.80 






3 


14.47 


2.1 


(0.80) 


^ Per 


M 


1 


6.76 


2.4 


0.75 






2 


8.67 


0.6 


(0.80) 






3 


10.26 


1.7 


0.62 






4 


15.51 


0.8 


0.65 


a Cam 


K 


1 


-4.61 


3.3 


1.44 






2 


0.18 


3.7 


1.60 


23 Ori 


K 


1 


21.72 


11.2 


0.89 






2 


24.42 


4.8 


[0.80] 


A Ori 


M 


1 


7.58 


1.2 


1.61 






2 


23.53 


1.3 


(0.80) 






3 


24.98 


1.7 


(0.90) 






4 


26.76 


4.2 


0.51 






5 


28.07 


1.7 


0.80 


t Ori 


K 


1 


8.84 


1.4 


[0.40] 


e Ori 


A 


1 


3.13 


1.3 


0.38 






2 


11.08 


0.9 


0.62 






3 


12.82 


0.9 


0.68 






4 


17.41 


0.2 


(0.50) 






5 


22.75 


0.3 


0.44 






6 


25.30 


1.0 


1.24 






7 


28.93 


0.3 


(0.50) 


C Ori 


M 


1 


-1.81 


1.7 


[0.50] 






2 


23.05 


0.8 


0.55 



Table 2 — Continued 



Star Obs^ Comp vq N 

(km s-l) (109 cm-2) (km s"!) 







3 


26.44 


0.9 


0.59 


K Ori 


K 


1 


16.76 


0.5 


(0.70) 






2 


20.94 


0.9 


(0.80) 


Ori 


K 


1 


7.64 


2.2 


0.86 






2 


10.73 


7.4 


(1.50) 






3 


15.36 


31.4 


2.23 






4 


19.47 


0.9 


(1.00) 


HD 72127 A 


A 


1 


13.13 


2.4 


0.60 






2 


16.64 


1.2 


0.53 


5 Sco 


A 


1 


-14.46 


0.4 


0.68 






2 


-11.24 


2.0 


0.90 


Sco 


A 


1 


-10.32 


1.4 


[0.80] 






2 


-8.92 


3.0 


0.72 


(T OCO 


A 

A 


1 


—6.00 


1.6 


0.63 






2 


-4.65 


1.5 


0.73 






3 


-2.65 


0.4 


0.25 






± 


— io.o4 


Q 1 


(^U.oU ) 






2 


-10.40 


4.4 


(0.80) 






3 


-8.00 


16.7 


1.15 


X Oph 


K 


1 


-11.21 


5.1 


0.67 






2 


-8.48 


2.1 


(0.80) 


COph 


M 


1 


-16.32 


0.5 


(0.80) 






2 


-14.65 


1.6 


0.69 






3 


-13.74 


0.8 


[0.50] 


HD154368 


K 


1 


-14.94 


6.5 


1.42 






2 


-5.03 


7.2 


(0.50) 






3 


-2.36 


6.4 


(0.50) 


PCyg 


K 


1 


-22.44 


3.9 


(1.00) 






2 


-19.84 


17.8 


0.84 






3 


-17.43 


3.8 


(1.00) 






4 


-10.35 


4.1 


(1.00) 






5 


-8.02 


3.7 


(1.00) 


a Cyg 


M 


1 


-20.86 


0.4 


0.84 






2 


-13.38 


1.3 


0.65 






3 


-9.94 


0.9 


0.28 






4 


-8.78 


0.4 


(0.70) 
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Table 2 — Continued 



Star Obs'' 


Comp 


Vq 


N 


fob 






(km s-l) 


(10^ cm-2) 


(km s-i) 




5 


-3.65 


0.5 


(0.60) 




6 


-2.69 


0.6 


0.33 


HD207198 K 


1 


-19.32 


3.4 


(1.00) 




2 


-15.72 


10.4 


1.66 




3 


-12.20 


8.3 


(1.00) 




4 


-8.92 


26.0 


1.16 


9 Cep K 


1 


-24.03 


1.8 


(0.70) 




2 


-18.03 


4.8 


2.05 




3 


-15.02 


13.0 


1.09 




4 


-10.79 


9.9 


1.59 


19 Cep K 


1 


-18.03 


6.4 


(1.50) 




2 


-14.33 


4.4 


(1.40) 




3 


-11.69 


7.5 


(1.40) 




4 


-8.22 


6.3 


[1.20] 


A Cep M 


1 


-34.67 


1.5 


(0.80) 




2 


-32.75 


1.1 


(0.60) 




3 


-30.71 


1.2 


(0.60) 




4 


-23.31 


0.8 


(0.80) 




5 


-17.56 


1.5 


(0.80) 




6 


-15.57 


2.0 


(0.60) 




7 


-14.06 


6.5 


0.63 




8 


-12.69 


2.7 


(0.80) 




9 


-10.80 


7.0 


0.89 




10 


-9.23 


3.3 


(0.80) 




11 


-5.79 


3.0 


1.25 




12 


-2.73 


7.3 


0.56 


1/ Cep M 


1 


-23.80 


2.4 


(0.80) 




2 


—22.32 


3.4 


(0.80) 




3 


-20.71 


6.9 


0.71 




4 


-19.14 


2.3 


(0.80) 




5 


-17.22 


12.1 


1.40 




6 


-14.38 


1.2 


(1.00) 




7 


-12.14 


2.1 


[0.40] 




8 


-11.40 


2.3 


[0.40] 




9 


-10.33 


34.2 


0.44 




10 


-9.10 


1.0 


(0.50) 


6 Cas K 


1 


-19.02 


7.3 


0.89 




2 


-10.85 


8.5 


1.24 



''Observatory: M = McDonald (FWHM = 0.53 km s'^); A = AAO 
(FWHM = 0.3 km s'^); K = KPNO (FWHM = 1.3-1.5 km s'^). 

''Values of b in square brackets were fixed to facilitate convergence 
of the fits, but are well-determined. Values of b in parenthese were also 
fixed, but are less well-determined. 
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Table 3. Abundances and Atomic Data 



Element Aq"- Xion T x lO"'' A(log ar x lO^^b T/or ag x 10"<i T/og x 10^ 
(eV) (s~^) (cm^s^-'-) (cm^'"*) (cm'^s"^) (cm~^) 



Li 


3.31 


5.39 


34. 




9.6 


35.4 






C 


8.55 


11.26 


20. 


-0.70 


8.3 


24. 


1.62 


12.3 


Na 


6.31 


5.14 


1.3 


-0.46 


5.8 


2.2 


0.93 


1.40 


Mg 


7.58 


7.65 


7.9 


-0.35 


5.9 


13.5 


1.18 


6.70 


Al 


6.48 


5.99 


390. 


-0.37 


15.4 


253. 






Si 


7.55 


8.15 


300. 


-0.39 


8.4 


357. 


1.16 


259. 


S 


7.27 


10.36 


78. 


-0.64 


10.5 


74. 


1.03 


75.7 


K 


5.13 


4.34 


5.6 


-0.49 


5.5 


10.1 


0.76 


7.37 


Ca« 


6.34 


6.11 


37. 


-0.42 


5.6 


66. 


0.85 


43.5 






11.87 


0.14 


-0.73 


27.9 


0.05 


2.34 


0.60 


Cr* 


5.68 


6.77 


80. 




(6.0) 


133. 






Mn 


5.53 


7.43 


14. 




5.45 


25.7 


0.75 


18.7 


Fe 


7.51 


7.90 


20. 


-0.41 


5.5 


37. 


0.75 


26.7 


Ni 


6.25 


7.64 


23. 




5.56 


41.4 







''Solar system abundances (logarithmic) from Anders & Grevesse 1989; Grevesse & Noels 1993 (H = 
12.0). 

''Pequignot & Aldrovandi 1986; photoionization rate T for WJl radiation field; radiative recombination 
coefRcient <Xr for T = 100 K. 

"^Difference in log F (edge minus center) for cloud of total Ay ~ 1.0, for grain model 2 of Roberge et 
al. 1981. 

''Grain recombination coefficient (per H atom) from Weingartner & Draine 2001; for T = 100 K., tj} = 
400 kV2 cm3. 

''Second entry is for Ca+. 

fp from Meyer & Roth 1990; Or estimated. 
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Table 4. Correlation Probabilities 



Y 


X 


Ntot^ 








Pc'^ 


Pk^ 


Ps" 


Ca I 


K I 


82 


19 


0.749 


0.000 




0.000 


0.000 


Fe I 


K I 


51 


18 


0.768 


0.000 


0.000 


0.000 


0.000 


S I 


K I 


19 


3 


0.926 


0.000 


0.000 


0.000 


0.000 


re 1 


T 

Ca 1 


oy 


23 


0.664 


0.000 




0.000 


0.000 


Ca I/Na I 


Ca II/H 


37 


13 


0.569 


0.009 


0.118 


0.036 


0.084 


Ca I/K I 


Ca II/H 


33 


7 


0.539 


0.012 


0.020 


0.007 


0.010 


Ca I/S I 


Ca II/H 


15 


1 


0.697 


0.006 


0.007 


0.043 


(0.086) 


Ca I/C I 


Ca II/H 


24 


5 


0.637 


0.004 


0.049 


0.114 


(0.123) 


C I/K I 


CH+/H 


26 


9 


0.273 


0.298 


0.016 


0.128 


(0.120) 


C I/S I 


CH+/H 


13 


2 


-0.132 


0.653 


0.779 


0.705 


(0.679) 


C I/Na I 


CH+/H 


33 


17 


0.197 


0.469 




0.034 


0.030 


CN/CH 


Ca I/CH 


68 


27 


-0.016 


0.582 




0.001 


0.001 


CN/CH 


Fe I/CH 


59 


32 


-0.045 


0.828 




0.114 


0.130 


CN/CH 


K I/CH 


57 


19 


-0.110 


0.524 




0.492 


0.571 


Ca I/K I 


CN/CH 


ryl 


22 


-0.081 


0.664 




0.000 


0.000 


Fe I/K I 


CN/CH 


45 


25 


-0.267 


0.263 




0.026 


0.021 



''Ntot = total number of points (detections plus limits); Nii^ = number of points 
with limits for one or both values. 

''Linear correlation coefficient and associated probability of no correlation, using only 
detected points (Bevington 1969); average values considering either X or Y as the 
dependent variable. 

"^Probability of no correlation from Cox's proportional hazard model, including limits 
(ASURV; Isobe et al. 1986). 

"^Probability of no correlation from generalized Kendall's tau method, including limits 

(ASURV; Isobe et al. 1986). 

"'Probability of no correlation from generalized Spearman's rank order method, in- 
cluding limits (ASURV; Isobe ct al. 1986). Values in parentheses are uncertain due to 
small number of points (<30). 
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Table 5. Correlation Fits 







log[Ar(Y)] 


= A + B X log[Af(X)] 








Y 


X 






RMS'^ 




Sample'^ 


Program^ 


Ca I 


K I 


3.11 ±0.68 


0.60 ± 0.06 


0.21 


62 


3/3 


R 






-1.66 ±0.32 


(1.00) 


0.23 


61 


3/3 


R 






3.71 ± 0.72 


0.55 ± 0.06 


0.20 


54 


3/1 


R 






-1.66 ±0.31 


(1.00) 


0.22 


53 


3/1 


R 


Fe I 


K I 


2.31 ± 1.09 


0.80 ± 0.09 


0.19 


32 


3/3 


R 






0.00 ±0.27 


(1.00) 


0.19 


32 


3/3 


R 






1.14 ± 1.58 


0.90 ± 0.13 


0.17 


26 


3/1 


R 






-0.10 ±0.25 


(1.00) 


0.17 


26 


3/1 


R 








0.66 ± 0.08 




01 


Of O 


A 


S I 


K I 


-4.40 ± 1.94 


1.52 ± 0.17 


0.16 


16 


3/3 


R 






1.62 ±0.24 


(1.00) 


0.17 


15 


3/3 


R 






-1.11 ± 1.82 


1.24 ± 0.16 


0.12 


12 


3/1 


R 






1.64 ±0.20 


(1.00) 


0.14 


12 


3/1 


R 






-2.68 ± 1.49 


1.37 ± 0.13 




19 


3/3 


A 


Fe I 


Ca I 


0.45 ± 1.47 


1.13 ± 0.14 


0.21 


36 


3/3 


R 






1.75 ± 0.29 


(1.00) 


0.21 


36 


3/3 


R 






1.03 ± 1.56 


1.06 ± 0.15 


0.18 


30 


3/1 


R 






1.66 ± 0.26 


(1.00) 


0.19 


30 


3/1 


R 


Ca I/Na I 


Ca II/H 


12.26 ±2.91 


1.75 ±0.40 


0.32 


24 


3/3 


R 






5.49 ± 0.49 


(1.00) 


0.35 


24 


3/3 


R 






n KQ -4- O OQ 


0.48 ± 0.25 




O 1 


of o 


A 


Ca I/K I 


Ca II/H 


11.45 ±2.52 


1.45 ± 0.32 


0.28 


26 


3/3 


R 






7.41 ± 0.40 


(1.00) 


0.28 


26 


3/3 


R 






3.56 ± 2.09 


0.58 ± 0.23 




33 


3/3 


A 


Ca I/S I 


Ca II/H 


16.22 ±4.32 


2.17 ± 0.57 


0.28 


14 


3/3 


R 






5.62 ± 0.46 


(1.00) 


0.33 


14 


3/3 


R 






4.81 ± 2.93 


0.90 ± 0.32 




15 


3/3 


A 




Ha TT /TT 
wa 11/ XI 


1 c 1 n -4- A on 


2.25 ± 0.55 


0.37 


1 Q 


Of O 


r> 
JTL 






3.92 ± 0.59 


(1.00) 


0.41 


19 


3/3 


R 






0.77 lb 2.41 


0.64 ± 0.27 




24 


3/3 


A 




CH+ /H 




0.19 ± 0.17 


0.30 


1 7 


o/ o 


Ti 
ix 






3.48 ± 0.31 


(0.00) 


0.32 


17 


3/3 


R 


C I/S I 


CH+/H 


2.09 ± 1.74 


0.02 ± 0.22 


0.26 


11 


3/3 


R 






1.90 ± 0.27 


(0.00) 


0.27 


11 


3/3 


R 


C I/Na I 


CH+/H 


2.07 ±0.96 


0.06 ± 0.12 


0.21 


15 


3/3 


R 






1.57 ± 0.21 


(0.00) 


0.22 


15 


3/3 


R 


CN/CH 


Ca I/CH 




< -2.5 


0.32 


41 


3/3 


R 


CN/CH 


Fe I/CH 




< -2.0 


0.32 


27 


3/3 


R 


CN/CH 


K I/CH 




< -16.0 


0.20 


36 


3/3 


R 


Ca I/K I 


CN/CH 


-1.93 ± 0.10 


-0.10 ±0.10 


0.20 


31 


3/3 


R 






-1.84 ±0.20 


(0.00) 


0.20 


31 


3/3 


R 


Fe I/K I 


CN/CH 


-0.34 ±0.13 


-0.27 ±0.14 


0.23 


20 


3/3 


R 






-0.11 ±0.14 


(0.00) 


0.23 


20 


3/3 


R 



''Uncertainties in A were determined in the fits when B was varied, but are the standard 
deviations of the set of values log{iV(Y)/[Ar(X)]B} when B was fixed. 

''Values in parentheses were fixed. 

'^RMS distance of points from best-fit line. 

"^Number of stars in sample. 

"n/m: n = 1 includes high resolution (FWHM < 2 km s~^) data only; n = 3 includes also 
lower resolution data; m = 1 excludes Sco-Oph and other "discrepant" stars; m = 3 includes 
all stars. 

'R = regrwt (no limits; weighted by I/ct'^); A = ASURV EM algorithm (limits; no weights). 
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Table 6. Relative Electron Densities^ 



Ion 


23 Ori SLV'' 


1^ Ori main'^ 


1 Sco main'* 


S Sco main° 


C Oph 


Li I 
Li II 

Tie 


< 9.70 
(11.4) 
<0.71/ •■■ 


(11.0) 
... / ... 


(11.8) 
... / ... 


8.92 
(11.8) 
0.047/ •■■ 


9.60 
(11.8) 
0.22/ ••• 


C I 
C II 

Tie 


14.98 
16.95 
0.258/0.214 


13.53 
(16.45) 
0.029/0.007 


13.91 
17.20 
0.062/0.004 


14.20 
17.38 
0.080/0.020 


15.51 
17.32 
0.373/0.077 


Na I 
Na II 

Tie 


13.35 

(14.4) 
0.200/0.093 


12.23 
(14.0) 
0.038/0.009 


12.50 

(14.8) 
0.056/0.010 


12.72 

(14.8) 
0.093/0.028 


13.85 

(14.8) 
0.251/0.023 


Mg I 

Mg II 

Tie 


13.80 
15.57 
0.227/0.206 


12.96 
15.24 
0.070/0.041 


12.98 
15.67 
0.137/0.025 


13.33 
15.60 
0.360/0.240 


13.80 
15.32 
0.404/0.239 


Al I 
Al II 
Tie 


< 11.53 
13.48 
< 2.8/ • • • 


< 10.7 
12.95 
< 1.4/ • • • 


... / ... 


< 10.55 
(13.2) 
<0.56/ ••• 


< 10.22 
(13.2) 
<0.26/ ■•• 


Si I 
Si II 

Tie 


11.79 

15.37 
0.094/0.081 


< 12.0 

15.14 
< 0.26/< 0.23 


11.40 

15.62 
0.108/0.028 


15.38 
... / ... 


15.34 
... / ... 


S I 
S II 

Tie 


12.81 
15.98 
0.050/0.036 


11.50 
(15.55) 
0.007/0.002 


11.6 
(16.4) 
0.006/0.001 


12.41 
(16.4) 
0.038/0.015 


13.81 
(16.4) 
0.191/0.059 


K I 

K II 

Tie 


11.50 
(13.15) 
0.228/0.096 


10.49 
(12.75) 
0.056/0.012 


10.88 
(13.55) 
0.109/0.015 


11.25 
(13.55) 
0.255/0.066 


11.85 
(13.55) 
0.203/0.016 


Ca I 
Ca II 

Tie 


10.20 
12.04 
0.955/0.955 


9.23 
11.57 
0.302/0.302 


< 9.7 
11.28 
< 8.7/< 8.7 


9.38 
11.50 
2.506/2.505 


9.46 
11.71 
0.372/0.370 


Cr I 
Cr II 

Tie 


12.64 
... / ... 


12.30 
... / ... 


12.98 
... / ... 


12.70 
... / ... 


9.8 
12.39 
0.34/ ••• 


Mn I 
Mn II 
Tie 


10.23 
13.08 
0.036/0.027 


< 10.0 
12.72 
< 0.05/< 0.03 


< 10.5 
13.33 
< 0.19/< 0.05 


... / ... 


< 11.3 
13.20 
< 0.32/< 0.21 


Fe I 
Fe II 

Tie 


11.34 
14.30 
0.040/0.040 


< 10.3 
13.98 
< 0.01/< 0.01 


11.07 
14.68 
0.045/0.045 


11.34 
14.51 
0.123/0.123 


11.23 
14.32 
0.030/0.029 


Ni I 


< 10.26 


< 11.7 


< 10.1 




< 11.4 
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Table 6 — Continued 



Ion 


23 Ori SLV^ 


^ Ori main'' 


1 Sco main*^ 


S Sco main« C Oph B* 


Ni II 


12.90 


12.44 


12.84 


12.91 




< 0.09/ • • • 


<7.5/ ■■• 


<0.08/ ■■• 


... / ... <i.3/ ... 



^Entries for ions are log N (cm~^). References for Li I, C I, Na I, K I may 
be found in WH; references for S I, Ca I, Ca II, and Fe I are in the appendix to 
this paper. Unless otlierwise noted below, values for other species are from fits to 
archival HST/GHRS spectra. Values in parentheses are estimated from N(Zn II) 
and typical depletions. Values tabulated here are for the main components in these 
lines of sight, and thus may be slightly different from the total line of sight values 
given in the references. Entries for Ue are cm~^, and are calculated via the T/ar 
values in Table 3 for T = 100 K and the average interstellar radiation field (except 
for 1 Sco and S Sco) . First rie is for radiative recombination only; second Ue includes 
grain-assisted recombination (see text for details). 

''23 Ori: SLV components at vq ~ 20 to 26 km s"! (Welty et al. 1999b). 

'^f Ori main: components at vq ^ 22 to 27 km s~^. 
'^1 Sco: components at vq ~ —9 to —5 km s""*^. 

"S Sco: components at vq ~ —14 to —11 km s~^. C II from Hobbs et al. 1982; 

AI I from Snow et al. 1986. 

Oph B: components at dq ~ — 16 to —12 km s~^. Some first ions adapted 
from Savage et al. 1992; C II from Cardelli et al. 1996; Al I from Shulman et al. 
1974; Or I from Meyer & Roth 1990. 
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Table 7. Pressure, Density, and Fractional Ionization^ 



Star log[Af(C I")] loglnuT]" nn log[Ar(Ca I)/Ar(Ca II)] 

Component log[Af(C I*)] fi^ (K) (cm-3) 

log[iV(C I**)] /2 



o Per'' 



15.59 
14.88 
14.16 



0.158 
0.030 



3.45±0.1 50 



55 



-2.08 



rie(Ca)<= ?ie(CNaK)f ne/nn^ 
(cm^"^) (cm~^) 



0.342 (0.218) 0.00396 
0.342 (0.089) 0.00162 



C Per'' 
main 



15.48 
14.25 
14.09 



0.054 
0.037 



2.85±0.1 55 



13 



-2.16 



0.304 (0.185) 0.01420 
0.304 (0.128) 0.00985 



23 Ori 
SLV 



14.93 
13.97 
13.15 



0.10 
0.015 



S.litO.l 100 



13 



-1.84 



0.955 0.229 0.01760 

0.955 0.134 0.01030 



A Ori'' 
main 



14.12 
13.74 
12.97 



0.28 
0.047 



3.65±0.2 45 



100 



-2.13 



0.284 
0.283 



(0.057) 
(0.005) 



0.00057 
0.00005 



C Ori 
main 



13.45 
12.66 
12.00 



0.13 
0.03 



3.4±0.2 (100) 



25 



-2.34 



0.302 
0.302 



0.041 
0.009 



0.00164 
0.00036 



v = 



13.16 
13.23 
12.84 



0.45 
0.18 



4.5±0.1 (100) 300 



-1.90 



0.832 
0.823 



0.162 
0.004 



0.00054 
0.00001 



1 Sco 
main 



13.50 
13.53 
13.21 



0.41 
0.20 



4.4±0.2 75 



335 



< -1.58 



<7.1 
<7.1 



0.062 
0.008 



0.00019 
0.00002 



5 Sco 
main 



14.15 
13.79 
13.08 



0.29 
0.056 



3.75±0.1 60 



95 



-2.12 



1.77 
1.77 



0.101 
0.028 



0.00106 
0.00029 



a Sco'' 
main 



14.03 
13.99 
13.40 



0.43 
0.11 



4.15±0.2 70 



200 



-2.18 



1.71 
1.71 



(0.067) 
(0.012) 



0.00034 
0.00006 



C Oph'' 
B 



15.33 
14.95 
14.38 



0.26 
0.074 



3.85±0.2 75 



100 



-2.25 



0.305 
0.304 



0.229 
0.033 



0.00229 
0.00033 



HD 185418'' 
main 



15.53 
14.45 
13.59 



0.081 
0.011 



2.8±0.2 100 



-2.32 



0.316 
0.316 



(0.161) 
(0.118) 



0.02680 
0.01970 



HD 192639'' 
main 



15.29 
14.53 
13.67 



0.134 
0.020 



3.2±0.1 100 



16 



(0.108) 
(0.048) 



0.00675 
0.00300 
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Table 7 — Continued 



Star 

Component 


log[JV(C 1°)] 
log[Ar(C I*)] 
log[iV(C I**)] 


h 


log[nHT]<= 


(K) 


(cm-3) 


log[Ar(Ca I)/iV(Ca II)] 


ne(Ca)« 
(cm-3) 


ne(CNaK)f 
(cm-3) 




A Cep'' 
main 




0.19 
0.07 


3.3±0.1 


70 


30 


-2.31 


0.254 
0.253 


(0.165) 
(0.049) 


0.00550 
0.00163 


v = -3 




0.30 
0.10 


3.8±0.2 


(100) 


60 


-1.54 


1.90 
1.90 






V = -35 




0.41 
0.36 


>4.8 


(100) 


>600 


-1.92 


0.793 
0.793 







^Unless otherwise noted, column density data are from same sources as in Table 6. 

''Practions of C I in first and second excited states. 
'^Estimated from graphs in Jenkins & Sliaya 1979. 

■^Prom H2 rotational excitation (Savage et al. 1977; Rachford et al. 2002); values in parentheses are assumed. 

'^Corrected for T < 100 K; WJl radiation field assumed (multiplied by 5 for Scorpius sightlines). First line for radiative recombination 
only; second line includes grain-assisted recombination. 

^Average value for C, Na, and K; corrected for T < 100 K; WJl radiation field assumed (multiplied by 5 for Scorpius sightlines). 
First line for radiative recombination only; second line includes grain-assisted recombination; values in parentheses are for entire line 

of sight. 

^Using rie for C, Na, K. First line for radiative recombination only; second line includes grain-assisted recombination. 

''C I fine-structure analysis from Wannier et al. 1999 (o Per, f Per); Zsargo 2000 (A Ori, a Sco); Lambert et al. 1994 ((^ Oph); 
Sonnentrucker et al. 2002 (HD 192639) and in prep. (HD 185418); and Jenkins & Tripp 2001 (A Cep). 
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Tablc 8. Fc I A3719 Component Stnictnrcs 



Star 


S/N 




Comp 




N 








(mA) 




(km s-l) 


(10^° cm-2) 


(km s-l) 


C Ori 


245 


<0.1 






<2.0 




HD 72127 A 


105 


<0.3 






<6.0 




p Leo 


105 


<0.3 






<6.0 




6» Car 


265 


<0.1 






<2.0 




<5 Sco 


285 


1.1±0.2 


1 


—13.42 


1.9 


(0 40") 








2 


-11.24 


19.8 


1.04 


/31 Sco 


345 


0.9±0.1 


1 


-10.36 


5.7 


1.10 








2 


—8.90 


12.1 


0.59 


<T Sco 


245 


2.0±0.2 


1 


-6.25 


22.7 


0.81 








2 


-4.43 


17.3 


0.88 


p Oph 


115 


11.9±0.4 


1 


-13.52 


1.9 


(0.50) 








2 


-10.33 


10.5 


0.32 








3 


-9.06 


64.2 


0.61 








4 


-7.87 


100.5 


[0.65] 








5 


-6.91 


98.2 


[0.65] 


C Oph*: 


175 


0.8±0.1 


1 


-14.75 


11.5 


0.58 








2 


-13.73 


5.0 


0.30 


^l Sgr 


125 


1.8±0.3 


1 


-9.64 


3.5 


(0.50) 








2 


-7.09 


7.3 


[0.80] 








3 


-5.53 


28.6 


0.73 








4 


-3.57 


2.3 


[0.30] 



''Uncertainties and limits are 2a. Chaffee 1974 found 3±1 mA for Per, 5±2 mA 
for o Per, and <1 mA for £, Per. Snow et al. 1986 found 1.0±0.8 mA for a Sco, <1 
mA for 5 Sco, and 13±2 mA for p Oph. Shulman et al. 1974 found 1.9 mA for ^ Oph. 

''Values of h in square brackets were fixed to facilitate convergence of the fits, but 
are well-determined. Values of b in parenthese were also fixed, but are less well- 
determined. 

'^Equivalent width and column densities from fit to profile shown by Barlow et al. 
1995, who gave b-values and velocities. 
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Table 9. Column Densities for Ca I, Ca II, Fe I, and S P 



tlD 


Name 


Wxl 


Ua i 


rtei 


oa ii 


net 


re i 


xler 


o 1 


rvei 


2905 


K Gas 


* 


9.92 


1 


12.70 


1 






13.38 


1 


5394 


'y Gas 


* 


<9.86 


2 


11.16 


2 










20041 






10.34 


3 






<12.56 


1 






21291 




* 


10.33 


1 






12.25 


2 






21389 






10.38 


4 






12.25 


2 






21483 






10.24 


4 






12.02 


2 






23180 


o Per 


* 


9.87 


1 


12.09 


2 


11.95 


1 


13.74 


1 


281159 


BD+31 643 




10.07 


3 














23466 






9.73 


6 














23630 


Tj Tau 




<9.51 


5 


11.19 


2 






<11.50 


2 


24263 


31 Tau 




10.15 


6 














24398 


( Per 


* 


9.76 


1 


11.96 


2 


11.78 


3 


13.40: 


1 


24534 


X Per 




9.95 


3 






11.82: 


1 






24912 


^ Per 


* 


9.72 


1 


12.27 


2 


<11.48 


1 


13.30 


2 


25940 


48 Per 




<9.46 


5 














27778 


62 Tau 


* 


<9.50 


2 


11.97 


1 


<11.60 


1 






29647 






10.05: 


3 






<12.00 


1 






30614 


CK Cam 


* 


9.85 


1 


>12.86 


3 


<11.78 


1 






34078 


AE Aur 


* 


10.32 


3 






12.23 


1 






34989 




* 


<10.03 


4 














35149 


23 Ori 


* 


10.20 


7 


12.10 


4 


11.34 


4 


12.81 


3 


36371 


X Aur 


* 


10.24 


3 






12.02 


1 






36822 


<j>^ Ori 


* 


9.97 


8 


12.42 


5 










36861 


A Ori 


* 


9.96 


1 


12.45 


1 










37022 


Ori C 


* 


<10.08 


5 


12.10 


2 






<12.40 


2 


37043 


t Ori 


* 


9.18: 


1 


11.93 


2 










37061 


NU Ori 


* 


<9.70 


4 






11.96 


1 






37128 


e Ori 


* 


9.67 


1 


12.04 


2 






<11.66 


2 


37643 






<9.56 


2 














37742 


C Ori 


* 


9.53 


1 


11.97 


2 


<10.48 


5 


12.10 


2 


37903 




* 


9.67 


9 


12.14 


6 


<11.60 


1 






38771 


K Ori 


* 


9.13: 


1 


11.94 


2 










41117 


Ori 


* 


10.62 


1 






11.90: 


1 






42087 


3 Gem 




10.45 


3 






12.05 


1 






43384 


9 Gem 




10.53 


3 














46056 






10.56 


3 














46202 






10.46 


3 






<11.90 


1 






46711 






10.64 


3 






<12.18 


1 






47839 


15 Mon 


* 


<9.56 


3 


12.33 


2 










48099 






10.28 


3 






<11.78 


1 






50064 






>10.88 


3 






12.40: 


1 






53367 






10.09 


4 






<11.90 


1 






54662 






10.14 


10 














62542 






<9.37 


11 






<11.66 


6 






63804 






>10.88 


10 














72127 


HD 72127 A 




9.56: 


1 


12.95 


7 


10.70 


7 






73882 




* 


<10.03 


10 


12.62 


1 


<12.62 


2 






80077 






10.46: 


10 














81188 


K Vel 




9.50 


12 


11.45 


8 











Table 9 — Continued 



HD Name WH*' Ca I<= ReF Ca II'' Ref^ Fe P ReP S I^ Ref^ 



89688 


23 Sex 


* 


9.84: 


6 


91316 


p Leo 


* 


<9.40 


5 


92693 






<10.03 


10 


93521 




* 


<9.86 


3 


100262 


o2 Cen 




10.03 


8 


110432 




* 


<10.03 


10 


116658 


a Vir 


* 


<8.86 


1 


141637 


1 Sco 


* 


<9.67 


4 


143275 


5 Sco 




9.38 


1 


144217 


/Ji Sco 


* 


9.65 


1 


144470 


LO^ Sco 


* 


9.34: 


8 


145502 


V Sco 


* 


9.92 


8 


146010 




* 


<9.86 


6 


147009 






<9.51 


9 


147165 


(T Sco 


* 


9.54 


1 


147343 






10.48: 


13 


147701 






<10.16 


13 


147888 


p Oph D 




10.33 


3 


147889 




* 


10.41 


14 


147933 


p Oph A 


* 


10.26 


1 


148184 


X Oph 


* 


9.86: 


1 


149038 


/I Nor 


* 


>10.94 


8 


149404 






10.46 


4 


149438 


T Sco 


* 


<9.16 


8 


149757 


COph 


* 


9.45 


1 


151804 




* 


<9.40 


8 


152236 


Sco 


* 


<9.16 


8 


154090 




* 


9.92: 


8 


154368 




* 


10.30 


1 


157056 


e Oph 




9.76: 


8 


157246 


7 Ara 


* 


<9.33 


8 


157841 




* 


10.07 


6 


159561 


a Oph 




<9.03 


5 


159975 


/I Oph 




9.56: 


5 


161056 






<10.16 


10 


164353 


67 Oph 


* 


9.81 


15 


165784 






10.05 


2 


166734 






>10.79 


3 


166937 


H Sgr 


* 


9.88 


8 


167263 


16 Sgr 




<9.76 


8 


167264 


15 Sgr 


* 


10.32 


8 


167971 






10.60 


2 


168076 






>10.92 


3 


169454 




* 


10.51 


10 


170740 






9.79 


3 


172028 






>10.75 


3 


179406 


20 Aql 




10.21 


16 


180554 


1 Vul 




<9.40 


5 


182568 


2 Cyg 




<9.26 


5 



11.54 


9 










12.03 


2 


<10.95 


5 






12.45 


2 










11.96 


10 










10.76 


2 










11.51 


2 


11.26 


7 


11.67: 


2 


11.62 


1 


11.34 


5 


12.41 


2 


11.73 


2 


11.26 


5 






11.67 


3 










11.95 


1 










11.40 


9 










11.86 


2 


11.60 


5 


12.63: 


4 






12.34: 


8 










12.34: 


8 










12.55 


1 










12.56 


9 






12.15 


1 


12.42 


5 


14.00 


5 


11.99 


1 






13.80 


5 


12.62 


1 










10.63 


2 










11.79 


2 


11.22 


10 


13.85 


6 


13.08 


1 










12.66 


1 


11.89 


11 


13.87 


7 


11.60 


1 










12.22 


11 










11.53 


2 










11.71 


3 












1 
1 


--'11 'TC 

<,ii. to 


1 
1 










12.15: 


1 










11.59 


5 










12.04: 


1 










12.37 


1 










12.14 


1 










11.62: 


1 










12.60 


1 










12.13 


1 


13.17 


8 


11.79 


3 
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Table 9 — Continued 



HD 


Name 


WH'' 


Ca P 


ReP 


Ca 11^ 


Refd 


Fe P 


ReP 


sif 


Reff 


183143 






10.53 


3 






<11.85 


1 






184915 


K Aql 


* 


<9.86 


2 














185418 






10.30 


3 


12.62 


1 


11.84 


1 


13.66 


9 


185859 






10.41 


4 














186745 






10.71 


2 














188209 




* 


<9.76 


2 














188220 






9.97 


6 














188350 






9.97 


6 














190603 






10.61 


2 














192639 






10.56 


4 






11.80 


1 


13.44 


9 


229059 


BD+40 4220 




10.51 
>10.76 


3 
3 






12.29 
<12.32 


1 






193237 


PCyg 


* 


10.52 


1 














197345 


a Cyg 


* 


9.61 


1 


12.29 


2 










198478 


55 Cyg 


* 


10.45 


17 


12.60 


3 


12.08 


1 






199579 






10.29 


3 






11.60: 


1 






200775 






<10.15 


18 














203938 






10.51 


4 






12.10 


1 






204827 






10.61 


3 






12.32 


1 






206135 






<10.00 


2 














206165 


9 Cep 


* 


10.47 


1 






11.97 


1 






206267 




* 


10.55 


3 






12.04 


1 






207198 




* 


10.68 


1 






12.22 


1 






207260 


ly Cep 


* 


10.83 


1 














207330 


TT^ Cyg 




9.40: 


5 


11.91 


3 










207538 






10.65 


4 














209409 


o Aqr 




9.56: 


5 


11.83 


3 










209975 


19 Cep 


* 


10.39 


1 














210072 






10.63 


4 






<12.62 


2 






210121 






10.16 


3 






<11.78 


1 






210839 


A Cep 


* 


10.58 


1 


12.76 


1 


11.81 


1 






212097 


32 Peg 




9.99 


5 


11.95 


3 










212571 


TT Aqr 


* 


9.86 


6 


12.10 


2 


11.20 


7 


12.08 


2 


217675 


o And 


* 


<9.16 


5 


11.73 


2 










218376 


1 Cas 

BD+63 1964 


* 


10.03 
>11.11 


3 
3 


>12.61 


3 


<11.85 
<12.00 


1 
1 






223385 


6 Cas 


* 


10.20: 


1 






<12.00 


1 







'^Entries are log( A'^) (cm ^ ) , and refer to the total value for the entire sightline. A colon indicates 
values for which the Icr uncertainty is greater than 25%. Upper limits are 3(7. 

''Sightlines included in the appendix to Welty & Hobbs 2001 are denoted by asterisks. 
"^Ca I: 1 = this paper; 2 = Allen 1994; 3 = Thorburn, priv. comm.; 4 = Welty unpubl.; 5 = 
White 1973; 6 = Penprase 1993; 7 = Welty et al. 1999b; 8 = Lambert & Danks 1986; 9 = Knauth 

et al. 2001; 10 = Grcdcl ct al. 1993; 11 = CardcUi ct al. 1990; 12 = Crawford 2002; 13 = CardcUi 
& Wallerstein 1986; 14 = Haffner & Meyer 1995; 15 = Federman & Hobbs 1983; 16 = Hanson et al. 
1992; 17 = Hobbs 1971; 18 = Federman et al. 1997. 

<^Ca II; 1 = Welty unpubl.; 2 = Welty ct al. 1996; 3 = Hobbs 1974, 1978: 4 = Welty ct al. 1999b; 
5 = Price et al. 2001: 6 = Knauth et al. 2001; 7 = Hobbs et al. 1991; 8 = Crawford 2002; 9 = 
Penprase 1993; 10 = Crawford 1995; 11 = Sahu et al. 1998. 

«Fe I: 1 = Thorburn, priv. comm.; 2 = Joseph et al. 1989; 3 = Chaffee 1974; 4 = Welty et al. 
1999b; 5 = this paper; 6 = CardeUi et al. 1990; 7 = Welty, unpubl.; 8 = Cardelli & Wallerstein 
1986; 9 = Haffner & Meyer 1995; 10 = Barlow et al. 1995; 11 = Snow et al. 1996. 

I: 1 = Welty et al. 1991; 2 = Welty, unpubl.; 3 = Welty et al. 1999b; 4 = Allen et al. 1990; 
5 = ZsaJTgo et al. 1997; 6 = Federman & Cardelli 1995; 7 = Snow et al. 1996; 8 = Hanson et al. 
1992; 9 = Sonnentrucker et al. 2002 and in prep. 



